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Chapter 1 

Introduction 



Quantum chromodynamics, or QCD, is the theory of the strong interaction. One of its most 
remarkable features can be seen in the way a hadronic system behaves, when it is either strongly 
heated or compressed. The system initially composed of individual baryons and mesons — 
particles whose substructure can at lower energies be studied only indirectly — transforms at an 
energy density of roughly 1 GeV/fm^ into a plasma of deconfined, though strongly interacting, 
quarks and gluons. This transition is today at the focus of considerable interest in the high energy 
physics community due to the fact that quark-gluon plasma, the phase of hot deconfined matter, 
is currently being produced in ultra-relativistic heavy-ion collisions at RHIC in Brookhaven. A 
similar large-scale experiment is also being prepared at LHC in CERN. 

The rapid progress in experimental heavy-ion physics has set an important challenge for 
theorists. One needs to have a solid understanding of the processes that take place in the nuclear 
collisions and in particular obtain accurate numerical predictions for the different quantities 
measured. Whether or not the plasma produced in the present-day experiments has had time 
to thermalize and reach an equilibrium state, it is clear that one of the most important and 
fundamental quantities describing the deconfined phase of QCD matter is the grand potential 
= —pV of quark-gluon plasma. Its value is of relevance both to the study of the evolution of 
the heavy-ion collision products in terms of ideal hydrodynamics, and to cosmology, as e.g. the 
cooling rate of the very early universe depends on the energy and entropy densities of its content. 

The purpose of the present paper is to study the perturbative evaluation of the QCD pressure 
at finite quark chemical potentials, i.e. keeping the net number densities of the different flavors 
arbitrary. It was originally written as the introductory part of a doctoral thesis regarding the 
subject and consisting of three articles ^ El EI 5 and we have therefore set as one of our main 
objectives to present a review of the theoretical and mathematical background behind this work. 
This includes following through the construction of the functional integral representation of the 
QCD partition function, deriving from it the finite-temperature Feynman rules of the theory, 
and finally in some detail outlining the calculation of the most recent 0{g^lng) result for the 
pressure at finite T and fi j2j. We will also analyze the renormalization scale dependence of the 
perturbative results and study their agreement with current lattice data, and in doing so will 
extend the treatment of the original articles somewhat. 

The paper is organized as follows. In chapter 2 the theoretical structure of QCD is reviewed, 
and some of the most fundamental results of finite-temperature field theory are presented. A 
special emphasis is given to the construction of the partition function of non-Abelian gauge 
field theories and to the derivation of generic rules for finite-temperature perturbation theory. 
Chapter 3 is then devoted to describing the history and present status of the perturbative 
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determination of the QCD pressure, and we analyze in particular how the description of the full 
four-dimensional theory can at high temperatures be reduced into studying three-dimensional 
effective ones. The special features encountered when trying to approach the limit of small 
temperatures and large chemical potentials in the perturbative expansion are highlighted as 
well. 

Chapter 4 contains a review and analysis of the results of Refs. ^Hj. We extend the discussion 
of Ref. 2j by studying the renormalization scale dependence of the perturbative expansion of 
the pressure and by making comparisons with the results of a very recent lattice computation of 
Allton et al. 7 . In chapter 5 we then finally draw conclusions, while the appendices contain an 
introduction to the analytic methods used in Ref. [2] to compute three-loop diagrams in finite- 
temperature QCD and in Ref. [S| to evaluate massive four-loop diagrams in three-dimensional 
field theories. Due to the highly technical and rather lengthy nature of this part we have wanted 
to keep it separate from the main body of the text. 

In the following we will always work in natural units, i.e. we set c = h = ks = ^- Other 
notational conventions will be fixed at the end of chapter 2. 
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Chapter 2 

Basics of QCD thermodynamics 



In this chapter we wih present a brief introduction to the general properties of QCD and to 
the basic principles of statistical mechanics. We will especially concentrate on deriving a proper 
functional integral representation for the QCD partition function, and end the chapter by listing 
the finite-temperature Feynman rules for the theory. More detailed presentations on the subject 
can be found e.g. in Refs. |51ini[Tn|. 



2.1 General properties of the theory 

Quantum chromodynamics is a renormalizable non-Abelian gauge field theory developed in the 
early 1970's to describe the observed features of the strong interaction ^21- Its symmetry 
group is the color group SU(3), but for the sake of completeness we will in this work consider 
a more general case of SU(A^) keeping the number of colors a free parameter. The basic 
constituents of the theory are then the N"^ — 1 gauge bosons, gluons, and the nj flavors of spin- 
1/2 fermions, quarks, whose bound states the hadrons are. The gluons belong to the adjoint 
representation of the gauge group and are denoted by where ^ = 0, ...,3 is a Lorentz and 
a = 1,...,A^ — 1 a color index. The quarks on the other hand transform according to the 
fundamental representation of SU(A^) and are here assembled into a single 4 Anj-component 
spinor ip with Dirac, color and flavor indices suppressed. In nature the number of flavors is of 
course six, corresponding to the up, down, strange, charm, top and bottom quarks. The last 
three have, however, masses above the GeV-scale, and as the particles with masses larger than 
the energy scale of interest (here a few hundred MeV) very efficiently decouple, we have in all 
practical applications nj <2>. In this work we will furthermore regard all relevant quark flavors 
as massless, making nj = 2 often the most realistic value of the parameter. 

One of the most fundamental properties of QCD is its asymptotic freedom |13l I14j . which 
states that at large energies the gauge coupling constant g of the theory approaches zero. This 
can be seen most easily from the running of the coupling as obtained from the leading order 
solution to the renormalization group equation, 

2 _ 247r2 

> - (lliV_2n/)ln(A/AQCD) ' 

where A is the renormalization scale and Aqcd ~ 150 MeV a free parameter corresponding 
to the characteristic energy scale of the theory. Asymptotic freedom implies that at very small 
distances the behavior of QCD tends to that of a free field theory making the use of perturbation 
theory feasible in the description of hard processes such as deep inelastic scattering. From 
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10 MeV -- 



color superconductivity 



308 MeV 



Figure 2.1: A schematic phase diagram of QCD jl5j in terms of the independent variables 
temperature and quark chemical potential ^ = ^-b/^- 

the thermodynamic point of view this means that at least in the limit of asymptotically high 
temperatures or chemical potentials one might expect a perturbative approach to be fruitful in 
the computation of the partition function. This is an important observation, and its validity 
will be examined in the following chapters. 

Before moving on to a more quantitative description of QCD thermodynamics let us briefly 
review the present understanding of the phase structure of the theory, which is summarized in 
Fig. 2.1. The (/u-T)-plane is divided into three parts corresponding to the three fundamental 
phases of strongly interacting matter. At temperatures below Tc ~ 160 MeV and chemical 
potentials below ~ 350 MeV, QCD matter appears in the form of hadrons in either a liquid 
or a gaseous form depending on the values of (j. and T. There is a first order liquid-gas transition 
line originating from the ground state of nuclear matter at (/i, T) = (308, 0) MeV and ending at 
a critical endpoint at T ~ 10 MeV. Above this temperature there is merely a smooth cross-over 
from the dilute gaseous phase into the more dense liquid phase. 

At temperatures and chemical potentials higher than Tc and //ci strongly interacting matter 
appears in a deconfined phase, where the quarks are liberated from their confinement. There 
is again a first order transition line dividing the hardonic and deconfined phases, which starts 
from approximately (350, 0) MeV and ends at a critical point at roughly (240, 160) MeV. If the 
matter is even more dilute than this, the transition is believed to become a cross-over. The 
details of this pattern are, however, strongly dependent on the number of flavors as well as the 
values of the quark masses (see e.g. Ref. jTH]). 

The deconfined phase of QCD is further composed of two distinct parts. At high temperatures 
strongly interacting matter forms a color-neutral plasma of almost freely moving quarks and glu- 
ons — the quark-gluon plasma — the equilibrium thermodynamics of which we will study in the 
present work. At low temperatures and large chemical potentials the properties of the system are 
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on the other hand completely different, as its behavior is dominated by an attractive interaction 
between two quarks lying on the Fermi surface. Analogously to the electromagnetic case these 
quarks form Cooper pairs, which now constitute the ground state of a third fundamental phase 
of strongly interacting matter, the color superconductor |17j . The superconducting phase has a 
rich internal structure on the (/i-T)-plane composed of e.g. a color-flavor locked and a crystalline 
phase, and it terminates in a melting transition at a temperature of the order of several MeV 
(for a review, see Ref. ^Hl)- The description of these phenomena is, however, outside the scope 
of our present treatment. 



2.2 Quantum statistical mechanics 

To prepare for the main topic of our presentation, the perturbative evaluation of the QCD 
pressure, let us first consider the thermodynamics of a general grand canonical ensemble, i.e. a 
system which can freely exchange energy and particles with its surroundings. Its most impor- 
tant equilibrium properties can be straightforwardly obtained as soon as one has computed the 
corresponding partition function, which is defined as the trace of the density operator taken over 
all the quantum states of the system, 

Z{V,T,{^ii}) = Tip = Trexp[-/3(ii-- J^^iiVi)]. (2.2) 

i 

Here H is the Hamiltonian, A', are the conserved number operators, /ij the corresponding chem- 
ical potentials and (5 = 1/T. Quantities such as the pressure, entropy and different particle 
numbers are available as first derivatives of Z through the fundamental relations 

P 

S 

Ni 

and the internal energy can furthermore be expressed as a linear combination of these. Different 
response functions, such as the susceptibilities, 

are as well directly obtainable as higher derivatives of the partition function, and thermal aver- 
ages of different operators can be determined through the formula 

{A) = ^Tr(M). (2.7) 

Apart from a few special cases, such as non-interacting systems, the partition function is 
unfortunately only rarely obtainable by analytically performing the trace of Eq. 1)2. 2() as an 
explicit sum over the different energy eigenstates. In field theories a useful general representation 
for Z can, however, be obtained by using an obvious analogy to the functional integral form of 
the transition amplitude in zero-temperature quantum field theory. The trace appearing in the 



ainZ 

dv ' 

d{T\nZ) 

Wf 

d\nZ 
T , 



(2.3) 
(2.4) 
(2.5) 
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definition of the partition function can namely be written as an integral^ over the states of the 
system as 



Z = I d4>{^\exp[-P{H-Y,f^iNi)]\(^), (2i 



and regarding the integrand as a transition amphtude over an imaginary time — i/3 we get for a 
neutral scalar field (r = it) 

P(/>exp[- /^dr / d^xC], (2.9) 



and for a charged Dirac field 



per. J 



ViljVi;exp[- /"^dr f d^x(£-^VV)]- (2-10) 

/ antiper. Jo J 

As in indicated here, the boundary conditions require the bosonic field to be periodic on the 
interval [0,/?], whereas the fermionic field is antiperiodic, 

0(r = O) = ^{t = P), (2.11) 
V'(t = 0) = -i;{T = P), (2.12) 
^(r = 0) = -i'ir = 13). (2.13) 

While the first of these relations is trivial, the origin of the latter two can be found from the 
Grassmannian nature of the integration variables in Eq. (j2.1()j) . 

Above we have adopted a notation, in which is regarded as a function of -0, 

V'^ = V^7o, (2.14) 

and all trivial normalization factors of the functional integrals are suppressed. The choice of 
the imaginary time formalism, on the other hand, implies that we from now on work exclu- 
sively in Euclidean metric with the identification xq —ixQ. This means that in all standard 
Minkowskian equations we must make the substitutions 

U 2 2 2 2 1 ir\ 

— xg Xj >■ Xg Xj — X^X^ I^Z.IOJ 

7i -in, (2.16) 

which renders all 7-matrices Hermitean according to the convention chosen. The time variable 
r will henceforth be always denoted by xq- 

It is well known that the analytic evaluation of functional integrals appearing in interacting 
field theories is in general impossible. Apart from a numerical approach, the only method 
available for us to compute the integrals of Eqs. 1)2.91) and 1)2. 1U() for a given Lagrangian is usually 
to resort to a perturbative approach, i.e. to expand the interaction part of the exponentiated 
action in a series in the coupling constant. This reduces the determination of the partition 
function to the evaluation of the connected vacuum diagrams of the theory with the main 
difference to the zero-temperature case being just the replacement of the momentum space loop 
integrals by sum-integrals — a direct consequence of the finiteness of the temporal direction 
in the coordinate space. In the next sections we will see an explicit example of the finite- 
temperature machinery in action, as we now move on to apply the above general discussion to 
the special case of QCD starting from the very basics, i.e. Eq. ()2.2)) . 

^To be specific, a functional integral over the space of normalizable and suitably smooth coordinate space 
functions. 
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2.3 Construction of the QCD partition function 



We define massless quantum chromodynamics by the bare Euclidean Lagrangian density 



Cqcd = \F;,F;, + i,Ip^l:, (2.17) 



where the field strength tensor and the covariant (fermionic) derivative read 

C ^ d^Al-d^Al + ar'^A^Al, (2.18) 
= d^-igA^ = d^-igA^T'^ (2.19) 

and where we have, as usual, defined 

IP = 7m^m- (2-20) 

The symbols T", a = 1, ...,N'^ — 1, denote the generators of the fundamental representation of 
SU(A^) and are N x N matrices in color space. They are used to define the usual antisymmetric 
and symmetric structure coefficients /'^'"^ and d"'"^ through the relations 

[r'^,r^] = ir^'^r^ (2.21) 
{r'^,r*} = id"''^r^ (2.22) 

and the antisymmetric coefficients furthermore act as the generators of the adjoint representation 
of the gauge group as 

{T^f" = if^^, (2.23) 
[r^r''] = if^^T^. (2.24) 

In this work we, however, stick to the convention of writing them out explicitly; e.g. the covari- 
ant derivative D will in the following always refer to that of the fundamental representation, 
Eq. (|2.19|) . Another convention we have chosen is to write all indices of the adjoint representation 
as superscripts. 

As can be straightforwardly verified, the QCD action is invariant under the local gauge trans- 
formations 

A^ = A^T'' ^ n-^Af,n + ^ {df,n-^) n, (2.25) 

il) n-'^ip, (2.26) 

where the transformation matrix J7 reads 

n = expligT^a""] (2.27) 

and a" are arbitrary smooth functions. This causes a problem when trying to evaluate the 
partition function of the theory, as a proper gauge condition must be imposed on the fields 
in order to avoid an overcounting of degrees of freedom: we must be able to choose a unique 
representative from all gauge orbits to contribute to the quantity. In the following we will briefly 
explain how this can be implemented in the functional integral representation of Zqcd starting 
from its definition 

Zqcd = TTexp[-P{H-Y,f^fNf)]. (2.28) 

/ 



7 



Here the sum runs over all relevant quark flavors, and the corresponding chemical potentials fXf 
may be assumed independent as long as we are dealing with energies sufficiently below the scale 
of the flavor-mixing weak interactions, mw- 

A convenient choice for the present purposes is to work in the Aq = (temporal) gauge, which 
leads to a relatively simple algebra but has the disadvantage of not fixing the gauge completely 
(see e.g. Refs. ^l^^). The Lagrangian now reads 

itemp = / ^ A" A" + I K.Kj + i^liD^i^ + V'^^ov} , (2.29) 

and defining the canonical momenta corresponding to the coordinates Ai in the usual way, 

= ^ = if, (2.30) 

we obtain as the Hamiltonian 

i^temp = j d'xH = j d''x\^Ii'iIi^i--^F^^F^^-ij-t,D,^-^^do^Y (2.31) 

Despite the gauge fixing, we still have the freedom to make gauge transformations independent 
of xq. This is a consequence of the fact that the Gauss law, 

G» = diFto + gr''A\F!^ + T'^^P^^I, = 0, (2.32) 

where is the generator of time-independent gauge transformations, doesn't appear in the 
Hamiltonian equations of motion 



= -TT-a^ M = (2.33) 



5H ■ _ 6H 

When evaluating the partition function we thus must impose Eq. (|2.32)1 explicitly on the states 
contributing to the result, which is most easily accomplished by including the operator 

P = [ VAexp[i(3 [ d^xA^G"] (2.34) 

JA(oo)=0 J 

into the definition of the partition function, where it acts as a projector onto the space of physical 
states. As is indicated above, only fields A, which vanish at spatial infinity, contribute to P. 
Following the treatment of Refs. |1U1 lllj . we now substitute P into Eq. (|2.28() in the form 

Zqcd = Jim^Tr rPexp[-/3(/7tcmp-5Z'"/^/)H) ' (^.35) 

which after the repeated use of the standard completeness relation and the renaming of the 
integration variable A as Aq leads to the functional integral representation 

Zqcb = I VA^VUiVijVi; (2.36) 

J Af_i per. 
ip antip. 



X exp 



^ dxo j d^x [n - ilildoA'i - iA^id^S'^^ - gr'''A^)U^ - (agT"^ + /x) v}] • 
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Here we have chosen periodic boundary conditions^ for the Aq field and have furthermore denoted 

fi = diag(/xi,/X2,...,/i„^). (2.37) 

The integral over the conjugate momenta in Eq. H2.36() is of a trivial Gaussian form and its 
subsequent evaluation leads us to the familiar expression 

Zqcd = / VA^Vi;ViPe^p[- dxo / d^x (£qcd - V'VV')] , (2-38) 

J A I, per. J 

ip anttp. 

analogous to the general results of Eqs. 1)2.9(1 and (|2.10() . This form of the partition function 
however still contains excess gauge freedom, as it — due to the periodicity choice made above 
— is invariant under SU(A^) gauge transformations periodic in the temporal direction, 

0(xo = (3) = ^{xq = 0). (2.39) 

This renders the functional integral singular due to the infinite volume of the corresponding 
local gauge group, ^ 

Vn = [ Vn = oo, (2.42) 

/ Q per. 

e su{]V) 

and implies that we must attempt to modify Eq. ()2.38|) in order to be able to factor out a piece 
proportional to this integral. 

To fix the gauge freedom remaining in the above form of Zqcd, let us now introduce the 
covariant gauge condition 

F-[A] ^ d^Al-r = 0, (2.43) 

where is an undetermined function of x. Following the standard Faddeev-Popov procedure 
we insert the functional 

A[^] = / V^5[F''[A^]], (2.44) 

/ f2 per. 

6 SU{iV) 

where AP" denotes the gauge transformed field of Eq. ()2.25|) . into Eq. (j2.38|) in the form AA^^ = 
1. Using the gauge invariance properties of the action 

S = j\xQ j A^x{CQCYi-^^li^), (2.45) 



■^This is based on mere computational convenience. It has been shown by Montonen (see Ref. |1T|) that in the 
temporal gauge the partition function is actually independent of the boundary conditions of A. 

■^Despite the compactness of SU(A^), the locality of the gauge invariance makes the volume of the group infinite. 
The integration measure is defined in such a way that in an infinitesimal region we may write 

VQ. = ]^dQa(x) (2.40) 

x,a 

and that the relation 

j VQ.f[Q\ = j VQ.f[Q.Q.'] (2.41) 
holds, where fi' denotes an independent gauge transformation (see Refs. |19ll2Ull^ l. 
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of the functional A and of the integration measure T^An, we then obtain 



Zqcd 



Afj, per. 
i/j antip. 



VAf^V^jViP exp[-S[A]] 



VAJ)il)V^^ A~\A] [ VfidlF^iA^jjexpl - S[A]\ 

i; antip. G SU{iV) 

[ VVL I ^ VA^Vi)Vi)/\-\A^]8[F''[A^]]eiLp\- S[A^W, (2.46) 

/ Clper. j A^^ per. ^ J 



G SU(iV) 



if) antip. 



where we in the final stage may clearly drop the O-integral as a trivial constant factor. This 
finally gives us the result 



QCD 



per. 
i/j antip 



VA^V^Vi) A"^ [A]5[F'' [A]] exp [ - S[A]\ , 



(2.47) 



from which all excess gauge freedom has been removed. 

What remains to be done is to express the functional in a more practical form, which 
can be achieved by writing it as a Faddeev-Popov determinant 



A'^[A] = det 



<5F'^(x) 



detM'^^ 



(2.48) 



F°=0 



a relation straightforwardly verifiable using Eq. (|2.4U|) . The determinant can, on the other hand, 
be expressed in the form of a standard Gaussian integral over Grassmannian 'ghost' fields r/, r^, 



det M' 



ab 



DrjDr] exp 



rj per. 



dxo I d'x I dyo I d'yt{x)M^\x,y)rj\y) 



,(2.49) 



where the integration variables obey periodic boundary conditions on the interval [0, (3] as a 
result of the invariance of the partition function under periodic gauge transformations. This 
appearance of anticommuting but periodic fields in the theory is thus merely a simple conse- 
quence of the choice we made for the boundary conditions of Aq. Other choices would have been 
possible, too, but would result in computational problems later. 

Inserting the above ghost field representation of A^^ into Eq. 1)2. 47(1 together with the explicit 
form of M"^ 



M''\x,y) = a^{(9^<5«'' + 5r^M^)5(x-y)}, 



(2.50) 



we have now successfully reduced the partition function into a well-defined form. It can be 
applied to practical calculations as soon as the remaining (5-functional has been removed, which 
is easily achieved by multiplying the functional integral by the factor 



exp 



1 

2? 



/3 



dxo / d\{nx)f 



(2.51) 



and subsequently integrating over the functions f"", on which the partition function clearly 
cannot depend. This leads us to the final result 



QCD 



per. 
antip 



VA^Vil^VipVfiV'q exp 



/3 



dxQ I d xCcQ 



(2.52) 
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where the effective Lagrangian density C^s is given by 

Ces = /:QCD + ^(5^^^)'-V'W + ^"(<9V' + ffr^M^5^)7?^ (2.53) 

From this point on this functional integral serves as the definition of the partition function for 
us. 

2.4 Feynman rules for finite-temperature QCD 

Not surprisingly, it is impossible to analytically evaluate the integrals in Eq. (|2.52|1 . The formula 
can, however, be used as the starting point of a perturbative treatment, which relies on the small 
value of the coupling constant at large energies and amounts to expanding the integrand in a 
power series in g. We will now demonstrate the procedure by first considering the case of a 
generic scalar field theory defined by its action S[(j)]. 

Let us start by dividing the action into a free quadratic part Sq and an interaction part Sj by 

S = So + Si. (2.54) 

An expansion of the partition function in powers of Sj, 

Z = fv(^e^p[-{So + Si)] = Y,^-^ l'v4>exp[-So]S^, (2.55) 

then reduces the determination of the functional integral into the evaluation of expectation 
values of different products of field operators. Using Wick's theorem, these functions can be 
represented in terms of ordinary coordinate space integrals, which contain in the integrands the 
two-point functions, or propagators, 

A{x-y) ^ I y"P(/>exp[-5o]0(x)0(y). (2.56) 

The integrals are graphically represented by the famous Feynman diagrams, in which each line 
represents the propagator A{x — y) and each vertex, which correspond to the different terms 
of the interaction part of the Lagrangian, a four-dimensional integration variable x. As can be 
straightforwardly verified (see e.g. Ref. 0), the computation of the partition function is equiva- 
lent to evaluating all connected vacuum diagrams that can be drawn using the propagators and 
vertices of the theory. In practice almost all diagrammatic calculations are however performed 
in momentum space, where the propagators represent momenta flowing through the graph. 
Returning to the special case of QCD, the separation of the action gives 

Si = 5eff-5o = j\xojd^x£eS-So (2.58) 

= 5^^dxo I d^xf^^r^'^r^'AlAlAl^Al + r'%d^Al)A'^Al - i^^M + P'^v'' A^d^r,'^, 

from where we can immediately read off the inverse propagators in coordinate space recalling 
the result of the Gaussian functional integral (as well as its fermionic counterpart) 

dxQ / d'^x0(x)A ^{x — y)(f>{x) 



1 

Z 



D(/)exp 



^{x)<l>{y) = A{x-y). (2.59) 
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The momentum space representations of the actual propagators can then be obtained by solving 
the corresponding Green's function equations, which is easily accomplished by means of Fourier- 
transforming. In the finite-temperature case a transform consists of a three-dimensional integral 
over the spatial part of the momentum and a sum over the discrete values of its zero-component. 
Suppressing the diagonal color parts of the propagators we easily^ get 

(A^W(P) = = (2.60) 



p2 I A'*^ ^' p2 

p2 



A^(P) = = -4 = -4' (2-61) 



\{P) = ^ = (2.62) 
where the zero-components of the bosonic (including ghosts) and fermionic momenta read 

(Po)bos. = 27rnr, n G Z, (2.63) 
(po)ferm. = (2n + l)^r - i^, n G Z, (2.64) 

respectively. The exact forms of these Matsubara frequencies are dictated by the periodicity 
conditions of the different fields in the temporal direction. 

The vertex functions of the theory are available from Eq. 1)2. 58p and read in momentum space 



a a 




-g^S^r'r^'iSaySps - S^sS^^) + ^ ^ ) ^ (2-65) 




Q 



igr'"'{Sa(i{P - Qh + ^f^iiQ - RU + Sa^iR - (2.66) 




gia{T%j, (2.67) 



igr'^Pa, (2.68) 



^Getting all the signs right in the Feynman rules can, however, be surprisingly tricky. See Ref. \'22\ for a 
thorough analysis of this point. 
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where the notation in Eq. (|2.(i5|) imphes a cychc permutation of the index pairs {b,(3), (0,7) 
and (d, (5). Apart from trivial factors of i, there are only two differences to the familiar zero- 
temperature Feynman rules, which both are due to the discrete nature of the zero-components of 
the momenta: the po (5-functions appearing in the vertices are replaced by Kronecker (5-symbols, 

<5(4)(p) ^ £si^)[p)5^^^,, (2.69) 

which ensures the conservation of momenta. In addition, conventional four-dimensional loop 
integrals are replaced by sum-integrals, 

po/IpoI 

where po stands for bosonic and {po} for fermionic momenta. Symmetry factors on the other 
hand remain unchanged at finite temperature, as does the rule of assigning a factor —1 to each 
fermion and ghost loop. This is again merely a trivial consequence of the anti-commuting nature 
of the Grassmannian variables. 

Having now defined a functional integral representation for the partition function and obtained 
the necessary Feynman rules, it is natural to start considering the computation of specific di- 
agrams. The analytic evaluation of sum-integrals is an interesting field of study in itself, and 
the techniques used in the calculations are constantly evolving. The traditional method relies 
on transforming the po sums into contour integrals by the Residue theorem, but this becomes 
increasingly difficult when the number of loops in a diagram increases. In the appendix B of 
the present paper we will introduce a more evolved scheme developed for the calculation of 
multi-fold sum-integrals at vanishing chemical potentials in Ref. [B] and generalized io n ^ 
in Refs. [^[15], and give specific examples of its application. The method consists of Fourier- 
transforming the spatial integrals into coordinate space, which makes it possible to evaluate the 
sums analytically and leads to the computation of one-dimensional hyperbolic integrals in the 
end. In the next three chapters we, however, restrict ourselves to only quoting the results for 
the different graphs and leave the explicit computations to the appendices. For now we instead 
concentrate on the physics involved in determining the partition function of QCD. 



2.5 Notation 

To finish the chapter we introduce some notation that will be in frequent use in the following. 
Using the generators of the fundamental representation of SU(A^) as well as the antisymmetric 
structure coefficients, we define four group theory factors (see e.g. Ref. [21] ) through 

^^^cd ^ jabcjabd ^ jygcd^ ^2.71) 

CpSi, ^ {T-T%, = ^^^<5.i, (2.72) 
Tpd"^ = TrT^T** = ^ (2.73) 

j^^cd _ ^abc^abd ^ 

N ^ ' 

where the trace is taken over both color and flavor indices. The dimensions of the adjoint and 
fermionic representations of the gauge group will henceforth be denoted by 

dA = d"-" = N'^-l, (2.75) 
dp = 5ii = dATp/Cp = Nuf. (2.76) 
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For certain frequently occurring combinations of special functions we will apply the abbrevi- 
ations 



C'(x,y) ^ 9,C(x,y), (2.77) 
^ c'(-n,z) + (-ir+k'(-n,z*), (2.78) 
K(z) = ^(z) + ^'(z*), (2.79) 

where C denotes the Riemann zeta function and ^ is the digamma function 

. (2.80) 

The properties of these functions will be examined in the appendix B of the present paper. 
The chemical potentials will henceforth usually appear in the dimensionless combinations 

fl = /i/(27rT), (2.81) 
z = 1/2 -ifi, (2.82) 

and in the context of considering the zero temperature partition function the following abbrevi- 
ation will be used: 

Y^fi} ^ /^'- (2.83) 
/ 

In sums over a single flavor index the subscript f in fif is usually suppressed. 

Even though we have so far been working strictly in four dimensions, we will throughout 
the following chapters apply dimensional regularization and use d = 4 — 2e. The momentum 
integration measure and the notation used for sum-integrals from here onwards read 

p/{p} 



T ^ j, (2.85) 

Po/{po} ^ 



where A = (47re ''')^/^A is the scale parameter of the MS renormalization scheme we have 
adopted and 7 is the Euler-Mascheroni constant. The excess factor of in Eq. (ITKll) wih be 
suppressed in the appendices. 
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Chapter 3 

Perturbative evaluation of the QCD 
pressure 

We are now ready to start applying the perturbative machinery to the evaluation of the QCD 
partition function. This is done assuming the energy density of the system to be large enough 
to guarantee that we work in the deconfined phase and that the coupling constant is sufficiently 
small for the expansion of the functional integral to be well-defined. In this chapter we review 
some of the most important work that has been performed during the last three decades to drive 
perturbation theory to new orders, but at the same time keep the main focus on the present 
status of the computations. The motivation of the multi-loop calculations will be addressed in a 
later chapter by quantitatively examining the convergence properties of the perturbative series 
of the pressure. 

3.1 A historical overview 

As soon as the status of QCD as the correct^ theory of the strong interactions had been es- 
tablished and the formalism necessary to deal with non-Abelian gauge field theories at finite 
temperatures developed, interest towards the perturbative determination of the QCD pressure 
began to rise. The leading order result for the quantity could be obtained immediately, as at 
vanishing coupling the theory reduces to a trivial one containing just non-interacting quarks and 
gluons. The pressure of such a system is given by the Stefan-Boltzmann law and reads 

Pqcd(9 = 0) ^ PsB = ^(A^'-l + ^)+-^r'E^/ + T^E/^/- (3-1) 

/ / 

The first non-trivial finite-temperature computation was on the other hand the determination 
of the order correction to this result, which was carried out independently by Shuryak and 
Chin in 1978 25 . This calculation consisted of evaluating the one- and two-loop vacuum graphs 
of the theory (see Fig. 3.1) and produced a sizable contribution to the pressure, 

^ / / ^ 

The T = part of this expression was, however, known already beforehand from a computation 
of Preedman and McLerran jlj, to which we will return soon. 

^According to our present understanding. 
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Figure 3.1: The one-, two- and three-loop Feynman diagrams contributing to the QCD partition 
function. Diagrams containing ghost hnes have been suppressed. 

Motivated by the magnitude of the 0{g^) term, it was natural to attempt evaluating the 
pressure further to three loops. It was however well understood, based on an analogy with 
scalar field theories, that conventional perturbation theory ceases to work at this order due to 
uncancelled infrared (IR) divergences originating from the ring diagrams of Fig. 3.1 c. The 
problem can be cured only through an all-orders summation of the IR sensitive graphs, which 
gives a finite contribution to the pressure, but at the same time produces the first signs of non- 
analytic behavior in to the perturbative expansion. The resulting order and ^^In^f terms 
for the high-temperature pressure were first obtained by Kapusta |22j and Toimela [JJj in 1979 
and 1983, respectively. 

The next 0{g^) and 0{g^) corrections to the pressure at vanishing chemical potentials were 
computed over a decade later by Arnold, Zhai and Kastening EH]) who made a significant 
contribution to the field of finite-temperature field theory by radically developing the compu- 
tational methods used to evaluate sum- integrals. This was soon followed by another important 
work of Braaten and Nieto j29j , who were the first to fully exploit the machinery of dimensional 
reduction |3n| I31j in this particular problem by showing how all terms in the perturbative se- 
ries non-analytic in g^ could be obtained using effective three-dimensional theories. Recently 
Kajantie, Laine, Rummukainen and Schroder jHSESl have used the same approach to drive the 
expansion of the pressure even one step further, order g^h\ g, in an impressive calculation involv- 
ing four-loop diagrammatics in the effective theories and an extensive use of the Form program 
p?5j designed for symbolic manipulation. This result has since then been generalized also to 
the case of finite chemical potentials by the present author in Ref. [2], where the perturbative 
orders g'^ through g^lng were simultaneously covered by generalizing the calculational methods 
of Ref. [31 to finite fi and taking full advantage of the existing 0{g^lng) result at /i = 0. 

As has been known for a long time [301) high-temperature perturbation theory runs into 
serious infrared problems again at order g^. This time they can be overcome only by applying 
genuinely non-perturbative methods, as can be most straightforwardly verified by investigating 
the IR structure of diagrams of the general type 




A simple power-counting (see e.g. Ref. [H]) reveals that all similar graphs with at least 5 loops 
contribute to the partition function starting at order , which necessitates the computation 
of an infinite number of diagrams in order to obtain the next term in the expansion of the 
pressure. Unlike in the case of the ring sum, no computational method has been found to enable 
the simultaneous treatment of the whole class of IR sensitive graphs, making it impossible to 
proceed further without resorting to lattice QCD. The effective theory approach has, however, 
helped one to isolate the non-perturbativity to a three-dimensional pure Yang-Mills theory, 
where lattice studies are considerably easier and less expensive than in the full four-dimensional 
QCD. We will return to the implications of this observation in the next section, where we present 
a brief introduction to the method of dimensional reduction. 

At zero temperature and large chemical potentials the history of perturbation theory is com- 
pletely different from the high-T case. The expansion of the pressure has been known to 0{g'^) 
already since 1977, when Freedman and McLerran !4 performed the summation of the ring 
diagrams in this limit. Since then there has been very little progress^ in the field, as the most 
interesting phenomena occurring in this region of the (/i-T)-plane have turned out to be of 
fundamentally non-perturbative origin [T7| . Due to this reason and certain computational diffi- 
culties, the region of large chemical potentials and small but non-zero temperatures has actually 
been left completely untouched since the early work of Shuryak and Chin [25j . Progress towards 
an order result for the pressure applicable throughout the deconfined phase of quantum 
chromodynamics is nevertheless currently under way 37 . 

Apart from the traditional perturbative approach, numerous other methods have been applied 
to the computation of the QCD pressure. These include most importantly four-dimensional 
lattice simulations (see e.g. Refs. [3l3HlinniEI3l^)) which have produced plenty of reliable results 
at small chemical potentials and temperatures close to Tc, and different resummation schemes 
of perturbation theory, such as the hard thermal and dense loop approximations |42 | I43 | H H I45j . 
While a detailed description of these methods is out of the scope of our treatment, we want to 
end the section by mentioning two computational schemes of special interest to us. 

Moore, Ipp and Rebhan have recently evaluated the QCD pressure non-perturbatively in 
the limit of a large number of flavors — where the theory simplifies significantly — both at 
vanishing chemical potentials [3^1 and at finite n [17j- In doing so they have been able to extract 
the coefficients of the perturbative expansion of the quantity up to order g^ |17j) which in the 
limit of large provides an accurate numerical check for the results of Refs. [21 132] and in 
addition gives an estimate for the magnitude of the yet undetermined 0{g^) contribution to 
the // = pressure. In Ref. jlH] the analysis of the perturbative results has been extended 
to cover the entire (/i-T)-plane, and one has in particular been able to quantitatively examine 
the breakdown of dimensional reduction in the 0{g^\ng) pressure 2^ in the limit of large /i/T. 
Remarkably, methods similar to those used in the large- A'^y context have also been used to study 
the non- Fermi-liquid behavior of the specific heats of the full theory at low temperatures |49j . 

Another very interesting approach to perturbative computations is the 'numerical stochastic 
perturbation theory' developed by Di Renzo et al. [HOI; which is based on finding numerical 
solutions to the Langevin equation that the quark and gluon fields are required to obey. The 
form of the equation determines the evolution of the fields in stochastic time and provides a 
natural way of expanding the solutions in power series in g. This method has recently provided 
important results that can be used to extract the non-perturbative 0{g^) contribution to the 
QCD pressure originating from an effective three-dimensional Yang-Mills theory. 

■^Recently, the result of Ref. |1| has been verified in an analytic computation [5], which removed all numerical 
uncertainties from the result. No new perturbative orders have, however, been obtained. 
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3.2 The pressure at high temperatures 



Let us now move on to take a closer look into the most recent perturbative computations of 
Refs. [2II22]. This part of the paper is divided into two sections corresponding to different regions 
on the (//-T)-plane. In this first section we cover the region of high temperatures and relatively 
small chemical potentials giving special emphasis to explaining how three-dimensional effective 
theories can be applied in the computations. In the latter one we then review the current status 
of perturbation theory in the limit of large chemical potentials and small temperatures, and 
outline some possible future developments. 

3.2.1 Dimensional reduction and effective theories 

If the temperature of a strongly interacting system is well above the critical temperature of 
the deconfinement transition and the chemical potentials are negligible in comparison^, the 
description of thermodynamics simplifies somewhat. This is due to a phenomenon called di- 
mensional reduction, which refers to the observation that the large distance static correlation 
functions of many field theories can at sufficiently high temperatures be obtained via effective 
lower-dimensional theories It is based on the reduction of the relevant degrees of freedom 
to the zeroth Matsubara modes of the bosonic fields, as we will demonstrate in the following. 
The method has played a central role in making the recent 0{g^liiig) perturbative computations 
of the QCD pressure possible. 

To introduce the concept, let us first consider a generic field theory, in which the masses 
of the different fields are of negligible magnitude in comparison with the temperature, and all 
chemical potentials are taken to be zero. Just as we have seen in the previous chapter, the 
functional integral representation of the partition function leads to the requirement of periodic 
and antiperiodic boundary conditions for the bosonic and fermionic fields on the interval xq G 
[0,/3]. This in turn causes their Fourier representations to be of the form 

00 

cPix) = T e™(/<„(x), (3.4) 

n=— 00 

where the bosonic and fermionic Matsubara frequencies read 

io^ = 2mTT, (3.5) 
a;f = (2n+l)7rr, (3.6) 

respectively. Assuming the kinetic part of the Lagrangian to be of the usual quadratic form, the 
momentum space propagators of the field modes are then always proportional to l/(a;^ +P^)) 
and it is obvious that the non-static modes gain effective masses that grow linearly with the 
temperature. It would therefore seem natural to expect some sort of decoupling to take place 
for these fields at high temperatures, which would leave the zero-modes (po of the bosonic fields 
as the true degrees of freedom. 

To verify that the decoupling indeed takes place, we refer to a well-known theorem of Ap- 
pelquist and Carazzone |53j . stating that for a renormalizable (zero-temperature) field theory 
containing two types of fields with masses mi <C m2, the Green's functions with typical momen- 
tum scales p <C m2 may up to corrections of the orders (mi/m2)", (p/m.2)", n > 1, be calculated 
using a Lagrangian with the heavy fields removed. In the new Lagrangian the coupling constants 

■^This implies fif <C 2nT for all flavors. For a more detailed analysis, see chapter 4 and Refs. ^1 151|. 
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have modified values and new, possibly non-renormalizable, interaction terms often need to be 
introduced in order to reproduce the results of the old theory. 

The application of the decoupling theorem to the computation of the finite-temperature par- 
tition function seems at first sight extremely straightforward, as we are dealing with diagrams 
with vanishing external momenta, and the light fields — which now clearly are the bosonic zero 
modes — are furthermore massless. Here we, however, must take into account also the possible 
generation of thermal masses for these fields, which leads to corrections to the partition function 
proportional to positive powers of the coupling constant. Assuming our theory to be asymp- 
totically free in the ultraviolet (UV), it is anyhow clear that with increasing temperature the 
static fields eventually become the dominant degrees of freedom, meaning that the non-static 
modes can be integrated out from the description of equilibrium thermodynamics at least in 
the framework of perturbation theory. This implies that we are left with an effective theory 
with a Lagrangian independent of xq, which makes it a genuinely three-dimensional one. In 
addition, the new theory generally has a simpler structure than the original one due to the lack 
of fermions, and practical calculations are often easier to perform there, as the tedious frequency 
sums disappear from the Feynman rules. On the downside, it has to be remembered that the 
effective theory parameters can usually be obtained only through the original theory. 

Returning to massless QCD, it is now obvious that the dominant degrees of freedom con- 
tributing to the high-T pressure are the different Lorentz components of the zero mode of the 
gluon field's Fourier expansion, 

oo 

= T exp[zu;^o]^^,n(x). (3.7) 

n=— oo 

They can up to a multiplicative factor be identified with the electrostatic scalar field Aq (x) 
and the magnetostatic gauge field Af (x) of a three-dimensional effective theory, electrostatic 
QCD (EQCD). The Lagrangian of the new theory is most easily obtained by first writing down 
the most general one respecting all the correct symmetries and then determining the different 
parameters through matching computations in full QCD. The eventual super-renormalizable 
Lagrangian density reads j51j 

Ce = ^ Tr 4 + Tr [A, ^o]' + m|Tr + ^ ^ /i^ Tr Aj] 

+ A«(TrA2)2 + Ag)TV4 (3.8) 

with 

i^^. = d^A'^-d^At + gEr'^A^A'^j, (3.9) 
A = di-igEAi, (3.10) 

and while there could in principle appear also higher-order, possibly non-renormalizable opera- 
tors in Eq. (|3.8j) . it can be easily seen that they would only contribute to the pressure at order 
g'^ or higher 32 . We have therefore neglected them in our analysis. 

The parameters appearing in the above Lagrangian are determined by demanding that the 
effective theory produces the correlation functions of the original one at distances 1 / (gT) and 
higher, which in the perturbative framework this leads to power series expansions in g'^ for 
them. The determination of the coefficients in the expansions, however, often becomes a highly 
laborious task in practise; see e.g. Ref. [1^1 for details. 
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Using the effective theory approach, we can now write the pressure of QCD in the form 

PQCD(T,/i) = pE(T,M) + ^lny"pA"^^Sexp{-SE}, (3.11) 

where the function is regarded as a parameter of the effective theory and is obtainable by 
computing the strict perturbation expansion'' of the pressure in the full theory. The pressure 
has here been divided into two parts corresponding to the contributions of different momentum 
scales: pE represents the hard scale 2ttT, which originates from the lowest non-zero Matsubara 
modes, while the partition function of the effective theory describes the effects of the soft scales, 
gT and lower, corresponding to the screening of the static gluon fields. All non-perturbativity 
naturally resides in the latter part. 

As soon as gauge fixing has been taken care of, the partition function of the effective theory 
can be determined by non-perturbative methods. One may, however, further simplify its compu- 
tation by noticing that the theory still contains two dynamical momentum scales, gT and g^T. 
The non-perturbativity manifesting itself at order g^ in the expansion of the pressure originates 
completely from the magnetostatic sector, which corresponds to the fields Ai and to the scale 
g^T proportional to their non-perturbative screening masses. This means that we can further 
divide the partition function of the effective theory into two parts by writing 

^ In^ VA^ VA^o exp {-Se} = PM + ylnJ VA^ exp { -5m} 

= PM+PG, (3.12) 

where the last non-perturbative term denotes the pressure of a new effective theory describing 
the magnetostatic sector (MQCD following the terminology of Ref. [22] )• This theory has been 
obtained by integrating out the Aq field and is defined by the Lagrangian density 

Cm = ^1^4, (3.13) 

Fij = diA';-djA'} + gur'''A'lA''j, (3.14) 

which is parametrized by the coupling constant gu ■ This parameter is determined in the electro- 
static effective theory (EQCD), and analogously to the case of pE, the function pM is obtainable 
through the determination of the strict perturbation expansion of the EQCD partition function. 
This is the last step in the series of effective theories, as all scales softer than g'^T are absent 
from MQCD due to its confining nature. 

3.2.2 The result to order 

In our present treatment we do not wish to go into great details in explaining how the parameters 
of the effective theories can be computed. Instead, we restrict ourselves to presenting the general 
structure of the results relevant for the computation of the QCD pressure, and start this by first 
inspecting the part originating from the hard momenta of order 2TrT. This scale enters the 
perturbative expansion exclusively through the parameters of EQCD, which can be shown to 
have the forms 1321 



PE{T,fi) 

rA-2^ 



,4 



^(«E 3 + 0(e) j +j^(Pei + 0{e) ) + 0{g') , (3.15) 



^Meaning a pure diagrammatic expansion, where no form of resummation has been apphed. It leads to a 
power series in g^. 
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2-loop 



3-loop 



4-loop 





Figure 3.2: The two-, three- and four-loop skeleton diagrams contributing to pu [SS]- The solid, 
wiggly and dotted lines stand respectively for the adjoint scalar Aq, the gauge boson Ai and the 
ghost fields. 
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(3.16) 
(3.17) 
(3.18) 
(3.19) 



where the a's and /3's are 'matching' coefficients with a priori unknown values. The a's are 
needed in determining the pressure up to order g^lng and have been computed in Refs. [51 1291 1^ 
1511 at /i = and in Refs. [51 at ;U 7^ 0. The most non-trivial one of them is the coefficient oes; 
which represents the contribution of three-loop diagrams to the strict perturbation expansion 
of the pressure (see Fig. 3.1 c). Its computation at arbitrary chemical potentials is explained in 
great detail in Ref. [2j and some of the techniques used in this calculation are also demonstrated 
in the appendix B of the present paper. The (3 coefficients, on the other hand, are not needed 
until order g^, and one of them, /3ei, still remains unknown. Not surprisingly, it will also be by 
far the hardest one to evaluate. 

The soft scale gT contributes to the pressure through the parameters of MQCD, which we 
write in the forms [51 132j 
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Figure 3.3: The one-, two-, three- and four-loop scalar master diagrams remaining after scalar- 
ization |33j . The solid lines correspond to massive and dotted to massless scalar propagators. 
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(3.22) 



As explained in the previous subsection, the above expression for pM has been obtained by 
computing the diagrammatic expansion of the EQCD partition function |291 I33j up to four 
loops, which was a long-standing and highly non-trivial problem. It was finished only very 
recently with the determination of the last two coefficients omi and /3mi by Kajantie et al. 
|33j . a task made possible by an efficient use ^1] of the Form program |35j designed for large- 
scale symbolic calculations. The general strategy of the computation was to first carry out all 
tensorial contractions in the diagrams and then reduce the number of necessary integrals into 
a small class of scalar master ones by applying an efficient algorithm developed by Laporta [S]. 
The analytic and numerical evaluation of Euclidean scalar integrals analogous to these is the 
topic of Ref. jSj, and it is also briefly dealt with in the appendix B of the present paper. The 
skeleton (two-particle irreducible) vacuum diagrams of EQCD are depicted in Fig. 3.2. and the 
scalar diagrams that remain to be computed after the scalarization process in Fig. 3.3. In the 
latter graphs the propagators have the usual scalar form l/{p'^ + m?'), where the mass parameter 
has the values m = uie and m = corresponding to the ^0 and Ai fields, respectively. The 
vertex functions in the graphs are simply constants. 

Finally, the ultrasoft scale g^T contributes to the QCD pressure through the partition function 
of the three-dimensional pure gauge theory MQCD. In terms of the parameters of the theory 
the perturbative expansion of this quantity reads 



Pg{T) 



TA 



-2e 



dACA 



+ o(/), 



(3.23) 



where the dependence on chemical potentials is of order g"^ or higher and where the magnetic 
mass iriQ is given by 



CAgl + 0{g^) 



(3.24) 



The coefficient oq appearing in Eq. (|3.23jl has been computed perturbatively in Ref. [HH, but 
Pg can only be determined using lattice methods due to the non-perturbative nature of pQ. 
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Its computation is a highly non-trivial task, as alone the conversion of the lattice results into 
continuum regularization requires complicated multi-loop calculations in lattice perturbation 
theory (see Ref. [32j ) . In this process the recent results of Ref. ^50^ may turn out to be very 
useful. 

Adding up the different parts of the full theory pressure and expanding in g, we are now ready 
to write down its perturbative result as 
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(3.25) 



Here g is the renormalized gauge coupling of the theory and the explicit 1/e terms cancel with the 
ones hidden in the matching coefficients. The result contains both an explicit and an implicit 
(through the running of g) dependence on the scale of dimensional regularization, A, which 
gets cancelled order-by-order but completely vanishes only through an all-orders summation of 
perturbation theory. This induces an arbitrariness to the result, which we may try to reduce by 
choosing the value of A in a way that would minimize the effects of the higher order contributions. 
This is, however, a highly non-trivial problem, and we will return to it in the next chapter. Here 
we merely note that even before determining the exact values of the coefficients /?ei and (5q, we 
can easily find the explicit form of the scale dependence of the order g^ term. This is achieved 
by applying the two-loop running of g, 
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to Eq. ()3.25|) and demanding that the A dependence of the result is beyond order g^ . The result 
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we hereby get for the A dependent part of the term reads 
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Combined with this result, the above expression for the pressure, Eq. ()3.25|) . represents the 
current status of high-temperature perturbation theory and in a clear way shows what is missing 
from the complete order result for the quantity. The matching coefficients a, which are 
necessary to obtain the pressure to order i^^ln g, are listed in appendix A, while two of the 
/?'s, /3ei and (5q, still remain unknown at present. For a more detailed derivation of the above 
results and the matching coefficients, we refer the reader to the original work performed in 

Refs. [2ii2niEaE3iniiin2]- 



3.3 The pressure at low temperatures 

A graphic analysis of the properties of the above high-temperature expansion of the QCD pres- 
sure will be carried out in the next chapter, where we will also quantitatively address the 
question, at which values of chemical potentials the dimensionally reduced result is applicable. 
It is, however, clear already beforehand that at least in the limit, where the temperatures are 
well below Tc but the chemical potentials large enough to keep us in the deconfined phase, 
the assumption of the temperature being the dominant energy scale can no longer be correct. 
This makes the reliability of the above result highly questionable, as it was derived assuming 
the zeroth Matsubara modes of the bosonic fields to be the dominant degrees of freedom. At 
low T this clearly is not the case, and the mechanism responsible for the screening of the IR 
divergencies is furthermore more involved than in the high-temperature case (see Refs. |Hini)- 
In the region of large ^i/T we therefore have no effective lower-dimensional theories to work 
with, which radically alters the computational methods available. We will now briefly describe 
the status of the perturbative calculations in this limit, but at the same time keep in mind that 
our treatment is restricted to the ordinary quark-gluon plasma phase. It cannot be applied to 
the study of color superconductivity, which in nature is the phase of strongly interacting matter 
near the T = line. 

At small but finite temperatures the only existing perturbative result for the pressure is of 
order g^ |25j . To determine the next 0{g'^lng) and 0{g^) terms in the expansion would require 
an explicit summation of the full theory ring diagrams^ to an arbitrarily high loop order, which 
leads to the computation of a very complicated logarithmic sum-integral. This has so far been 
achieved only in the limits of large T [22] and T = but a project aimed at determining 
the 0{g^) pressure at all values of fi/T is currently being completed 37 . As has been argued 
in Ref. [2], the quantitative difference between the outcome of this new computation and the 
present dimensionally reduced result may, however, turn out to be small in most parts of the 
(li-T) -plane. The motivation for this claim will be studied in the next chapter. 

At exactly zero temperature several simplifications occur in the computation of the pressure 
[Hl^: sum-integrals reduce to ordinary four-dimensional integrals and, more importantly, only 
diagrams containing fermion loops give non-zero contributions to the ring sum. This is due to 

^Those containing the one-loop giuon polarization tensor (see Fig. 3.4 d) as loop insertions. 
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d) n^i.(P) = -w^^^^ H — i'^ j — \ — — h 

Figure 3.4: a) The fermionic part of the one-loop gluon polarization tensor divided into its 
vacuum (T = fi = 0) and matter (vacuum-subtracted) parts. 

b) An IR safe diagram contributing to the pressure at T = 0. 

c) The generic form of the zero-temperature ring diagrams. 

d) The whole one-loop gluon polarization tensor. 



the fact that apart from the scale parameter, the only dimensionful numbers appearing in the 
diagrams at T = are the chemical potentials, which enter exclusively through the fermionic 
propagators. When the temperature is zero, the scale-free gluon (are ghost) diagrams simply 
vanish under dimensional regularization. 

Following the treatment of Ref. the zero-temperature pressure can up to order g'^ be 
observed to receive contributions from three classes of diagrams. These include the infrared safe 
fermionic one-, two- and three-loop graphs of Fig. 3.1, the single IR safe diagram of Fig. 3.4 b. as 
well as the ring sum of Fig. 3.4 c, where the grey circle denotes the vacuum-subtracted part of the 
fermionic one-loop gluon polarization tensor. The values obtained for the IR safe diagrams at 
arbitrary T and in Ref. can be straightforwardly continued to the limit of zero temperature, 
while the ring sum must be separately computed at T = |21 E] • Adding the different parts 
together, we then obtain as the final result of the T = pressure in the MS scheme 
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where the function F is defined by 



2;.^ 5: In 4 + I 5: I (/., - ^,)^ln 
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and the constant 6 has the approximative value 5 ~ —0.85638320933... This constant possesses 
a simple one-dimensional integral respresentation 0, the analytic value of which we however 
have not been able to obtain even using the PSLQ algorithm 
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Chapter 4 

Analysis of the perturbative results 



The correctness of most of the terms in the 0{g^lng) perturbative result obtained for the QCD 
pressure in the previous section has been verified in several independent computations in the 
limits of n = T = |^ and large Uf The practical value of the result is, however, 

solely determined by its predictive power, i.e. in which parts of the (//-T)-plane we can trust 
perturbation theory and, more specifically, dimensional reduction to be applicable. As there is 
no experimental data available for the pressure, this question is a very difficult one to address. 
There are nevertheless several methods we may try to apply. 

Four-dimensional lattice simulations (see e.g. Refs. 131^1231101) have during the past decade 
produced a number of accurate results for the pressure, but are even today restricted to small 
chemical potentials, ^ T, and to temperatures relatively close to the critical one, T < 5Tc. 
This is due to the fact that at finite n the exponentiated action appearing in Eq. (|2.52|) becomes 
complex and thus spoils importance sampling, while at large values of T the lattice treatment 
rapidly becomes more and more time-consuming. The lattice computations may nevertheless 
provide a good check of our results in the region of relatively small T and fi, but it is a priori 
not at all clear whether perturbation theory can be meaningfully applied so close to the phase 
transition line. The perturbative expansions have after all been derived assuming the coupling 
constant to have an in principle arbitrarily small value, which is only guaranteed to be the case 
in the limit of asymptotically large temperatures or densities. 

Another useful way to investigate the reliability of perturbation theory — this time on the 
whole (/x-T)-plane — is to analyze the convergence properties of the perturbative expansions 
themselves. Assuming that a consistent all-orders summation of perturbation theory would lead 
to a correct result for the pressure, the question of determining the accuracy of an expansion 
carried out to a specific order in g can naturally be reduced to analyzing the rate at which 
the corresponding series converges. This can, on the other hand, be often easily and efficiently 
examined through only the first few terms of the expansion, as we will soon observe. 

In our analysis of the perturbative results we will from now on restrict the treatment to the 
physical case of = 3. We find that the most illustrative way to proceed is to plot the results as 
functions of temperature and chemical potentials and then to analyze the graphs in a qualitative 
fashion. This is done separately for three different cases: first for the pressure at zero chemical 
potential, then for the //-dependent part of the pressure at high temperatures, and lastly for the 
zero-temperature pressure. In each region we investigate the convergence properties and scale 
parameter dependence of the results and for the first two we also provide a comparison with 
recent lattice results. The analysis of the /i-dependence of the high-temperature pressure is of 
special interest to us, as there we have both the pressure and the quark number susceptibilities 



26 




Figure 4.1: Left: the order cf" perturbative result for the nj = pressure plotted for different 
values of the 0{g^) coefficient. Right: the result plotted to different perturbative orders. The 
graphs are from Ref. 32 and the lattice data from Ref. j59j . 



(see Ref. [B) to compare with the lattice data. 

When investigating the plots of the following sections, it should be noted that while in Sees. 4.1 
and 4.3 we consistently apply the one-loop result, Eq. 1)2. 1|) . for the gauge coupling constant g, 
this will not be the case in Sec. 4.2. There we will in analogy with Refs. ^ [2] use a three- 
loop result for the quantity, as we have wanted to make the comparison with the figures of 
these papers as straightforward as possible. The effect of this choice on the plots is numerically 
relatively small, but the question is nevertheless in principle important. Using a three-loop result 
for g is somewhat inconsistent in light of the fact that one only needs the leading order term 
of its running to verify that the renormalization scale dependence of the perturbative result, 
Eq. (|3.25l) . cancels at order g^\ng. 



4.1 The pressure at vanishing density 

At zero chemical potential the order g^\n.g expression for the QCD pressure was first derived 
in Ref. |S1]. The analysis of this result is complicated by the undetermined 0{g^) term in 
the expansion, which induces an arbitrariness to the numerical factors that appear inside the 
logarithms of the g^lvig part. This has a considerable effect on the behavior of the result, 
and until the full 0{g^) contribution has been computed, the best we can do is to treat the 
coefficient of this term as a free parameter. This has been done in Fig. 4.1 where the 
pressure has been plotted as a function of temperature at n/ = using the value^ Aopt = 6.742T 
for the scale parameter. It is seen that by optimizing the value of the unknown g^ part, the 
perturbative results can be made to cleanly approach the lattice datapoints even at temperatures 
surprisingly close to Tc, but that the effect of varying the coefficient of the 0{g^) term is so 
large that the predictive power of perturbation theory remains regrettably weak at present. 
As has been conjectured in Ref. ^^j) there is nevertheless reason to believe that once the full 

^This value has been obtained by optimizing the convergence of the one-loop (NLO) expansion of g^. See 
Ref. E?! for details. 
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Figure 4.2: Left: the different terms of the perturbative expansion of pM + Pg at Uf = 
normahzed by Tg^ = T'^{g^ + 0{g^)) [SI]. Right: the scale dependence of the = pressure 
(the g^{hi{l/g) + 0.7) case) at n/ = Lattice data is again from Ref. [5^] . 



order g^ computation has been finished, the situation will improve significantly. At this order 
the magnetostatic sector namely enters the expansion of the pressure non-perturbatively, after 
which one has accounted for the dynamics of all the relevant energy scales — 2itT, gT and g^T. 
The convergence properties of the ;U = expansion have been analyzed in detail in Refs. [291 
[Sni with the result that the slowest convergence has been traced back to the contributions 
of the soft momentum scales entering through the effective theories. This is due to the simple 
fact that while for hard momenta of order 27rT the expansion parameter is ^^/vr^, the situation 
is completely different for the soft scales: for the color-electric modes the expansion is in powers 
of g/i: and for color-magnetic modes there is no perturbative series at all. This part of the 
pressure, or more specifically pM +pg-, is studied in Fig. 4.2 a, where the absolute values of the 
different terms of the function are plotted appropriately normalized. The general tendency in 
the plot is clearly favorable in the sense that at sufficiently large energies the magnitudes of the 
different terms keep decreasing as we move on to higher orders. This, however, only happens 
at very high temperatures, and the convergence of the perturbative expansion close to Tc is far 
more questionable. 

The scale dependence of Eq. (|.S.25|) is an equally important issue to address, as it induces an 
ambiguity to the results, which in principle is at least as serious as that originating from the 
unknown order g^ term. There are, however, several methods we may try to use to reduce this 
effect, of which the most widely used is due to Brodsky, Lepage and Mackenzie [H^. The BLM 
criterion states that the scale should be chosen in such a way that in the highest perturbative 
order it makes the term proportional to the highest power of n/ vanish. The scheme has 
been applied to the expansions of the parameters of the effective theories in Ref. [5^ with the 
result that the optimal scale has been observed to vary roughly between vrT and ^ttT. This 
seems physically very reasonable in light of the fact that the lowest non-zero bosonic Matsubara 
mode 2'kT falls within this interval and that with this choice the magnitudes of the logarithms 
appearing inside the matching coefficients remain small. 

In Fig. 4.2 b the renormalization scale dependence of the pure gauge pressure (at /x = 0) 
is analyzed by plotting the function at A = Aopt/2 and A = 2Aopt. It is clearly seen that 
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Figure 4.3: Left: AP|„^.=o plotted for different values of fi. Right: the different perturbative 
orders of AP{fi = 0.44Tc)|„y=o normalized by the free theory result. The lattice data is from 
Ref. 41 , and the result Tc/A(^|„,=o = 1-15 from Ref. |n2] has been applied. 



while at high temperatures the effect of the scale variation rapidly decreases, it does produce 
a sizable ambiguity to the results at temperatures close to Tc. This implies that even after the 
determination of the order contribution to the pressure, the perturbative result will always 
contain one free parameter that can in principle be adjusted to match the lattice data at low 
temperatures. There have been attempts to improve the situation through different resummation 
schemes, where the contributions of the soft scales have been treated without truncations j56j . 
but a detailed analysis of this work is outside the scope of our present treatment. 



4.2 Effects of small but finite chemical potentials 

Moving on to study the chemical potential dependence of the high-temperature pressure [2], 
we have a large selection of lattice data available for comparison. The pressure itself has been 
computed at small values of fj, in Refs. [3 I41j . and at ^ = its higher derivatives with 
respect to the chemical potentials, i.e. the quark number susceptibilities, are available from 
Refs. [3 1^ IMl lEOj- The compatibility of the perturbative results and the mainly quenched 
lattice data of Gupta et al. |^ EOl ESI has been analyzed in detail in the papers ^ 12 , but here 
we will extend the analysis by comparing the behavior of the perturbative expansion also with 
the Uf = 2 work of Allton et al. [7] . 

The quantities we analyze in this section are the deviation of the pressure from its /x = 
value, 

AP(r,^) = Pqcd(T,^)-pqcd(T,0), (4.1) 
as well as the quark number susceptibilities 

In Ref. |41j these functions were computed on the lattice in the quenched approximation, where 
no sea quarks are present. Prom the perturbative expression of the pressure, Eq. ()3.25|1 . the 
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Figure 4.4: Left: XuulT'^\nj=Q plotted for different values of A corresponding to the expected 
effect of the 0{g^^ term. Right: the expansion of the quantity at A = to different perturbative 
orders. The lattice data is from Ref. p^ . 
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Figure 4.5: Left: the scale dependence of AP|„^=o at /.t = 0.44Tc. Right: the scale dependence 
of XuulT^\nj=^- The lattice results are from Refs. [ITl IH^ . 



corresponding results can on the other hand be obtained by first computing the susceptibilities 
at arbitrary nj, then taking the formal limit nj = and in the end integrating the result 
back to obtain AP (see the discussion in Ref. We have compared these results for AP 

and XtiU) the linear (second-order) diagonal susceptibility, in Figs. 4.3 and 4.4, which show 
remarkable agreement even at small values of T . In addition to this, we observe a considerable 
improvement in the convergence properties of the perturbative expansions in comparison with 
the zero chemical potential case. This time also the scale dependence of the results, shown in 
Fig. 4.5, remains almost entirely within the error bars of the lattice data already at T = ZTc- In 
studying the above figures it should in particular be noted that unlike in the = case we have 
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Figure 4.6: Left: Xuu/T'^\nf=2 plotted to different orders together with lattice data [HO]. Right: 
lattice results for different susceptibilities |7j (for definitions, see the original paper). In the first 
plot A = 8.112T 57 and T^/ hr^\nf=2 = 0.49 [HI] have been used; in the second Tq = Tc- 




Figure 4.7: Left: AP/T^|„j.=2 plotted for different ratios of ^/T", i.e. not for fixed ix/Tc- Right: 
lattice results for the same quantity from Ref. 1 7^ . 



here in no way optimized the coefficient of the g term in the perturbative result, but instead 
simply used the 0{g^hig) part of Eq. (|.S.25|) . 

Figs. 4.3 to 4.5 reveal also another very interesting feature of the quantities AP and Xuu'- 
when normalized by their respective free theory values, the corresponding curves practically fall 
on top of each other. This observation, which at first instant seems curious, is explained by 
the fact that at small values of /i the behavior of the quantity AP is almost solely determined 
by the leading order term of its Taylor expansion in /x^, the linear diagonal susceptibility. The 
linear non-diagonal susceptibility Xud is of negligible magnitude in comparison with Xuu, and 
the contribution of the non-linear susceptibilities to AP is suppressed by a factor , which 
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Figure 4.8: p/pss plotted on the (/i-T)-plane for = 2. 

for = 0.44Tc, T > Tc is very small. In the present work we will therefore not consider the 
non-linear and non-diagonal susceptibilities at all, but refer the reader to Refs. I52j. 

The physically more interesting case of two quark flavors with = l^d = has been considered 
on the lattice in Refs. [THHO]. For the linear diagonal susceptibilities the compatibility of these 
results and their perturbative counterparts is analyzed in Fig. 4.6, which shows good agreement 
at least with the data of Ref . 60 . Taking into account the relation 

Xuu = ^ (Xq + 4xi) (4.3) 

between the definitions of Refs. |2] and j7], the outcome of Ref. can, however, also be seen to 
be fit in very nicely with the perturbative result, with the exception of temperatures very close 
to Tc- A similar behavior is observed for the pressure in Fig. 4.7, where the quantity AP is 
plotted for different fixed values of /x/T. The good agreement of our predictions and the lattice 
results of Ref. j2j, obtained using an unphysically high value 0.4r for the fermion masses, implies 
that the dependence of the results on quark masses should be of negligible magnitude. 

4.3 The zero-temperature limit 

While we have just witnessed the dimensionally reduced perturbative expansion of the pressure 
do exceptionally well in all lattice and convergence tests especially at finite chemical potentials. 
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Figure 4.9: The 0{g'^) result for p/psB\nf=2 plotted for temperatures below Tc. 

it is — as we pointed out already in chapter 3 — a priori not known how far away from the 
line ^ = we may continue to use it. This question can again be ultimately answered only 
through comparison with experimental data or with the results of a genuinely non-perturbative 
first-principles calculation. In the absence of both alternatives we must here resort to other 
means and choose to follow and extend the treatment of Ref . . 

Let us first consider a strongly interacting system at a temperature well above the critical one, 
which ensures that we are working in the deconfined phase. As has been argued in Ref. |51j . 
we expect dimensional reduction and hence Eq. (|3.25|) to be applicable at least as long as the 
chemical potentials satisfy /xj < 4T for all flavors. There is, however, no specific reason why we 
should expect the method to break down even at this point^, and it is therefore worthwhile to 
investigate the behavior of our result also at relatively high /u/T. This has been done in Fig. 4.8, 
where p/pss has been plotted as a function of /z and T in the case nj = 2, fiu = ^J-d = fJ- with 
the scale parameter (somewhat arbitrarily) reading 

A = 2^TT^/TT2p?. (4.4) 

We observe that the quantity stays below its free theory result 1 everywhere and only slowly 
approaches it when the values of the temperature and the chemical potential increase. In the 

■^This will happen inevitably only when the chemical potentials reach values /i ~ 2-kT / g making the three- 
dimensional mass parameter to be of the same order of magnitude as the hard scales. 
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Figure 4.10: The behavior of p/psB.\nf=2 at 0{g'^) near the T = hne. 

region of low T and ji the result on the other hand rapidly decreases reaching a value close to 
zero when T = Tc and = 0. There the applicability of perturbation theory is, however, anyway 
highly questionable due to the large value of the coupling constant. 

Moving down along the temperature axis, it is clear that the dimensionally reduced result 
ceases to be applicable at some point. If we let the chemical potentials to be high enough to keep 
us in the deconfined phase, we may, however, keep inspecting the behavior of the function and 
especially compare its order g'^ part with the zero-temperature result of Eq. 1)3.27(1 . Continuing 
to work at nj = 2, we observe from Fig. 4.9, where the pressure has been plotted to 0{g'^) at 
temperatures below Tc, that the perturbative result slowly decreases as a function of T all the 
way down to T/ (j, ~ 1/40. At even lower temperatures than this the unphysical nature of the 
result^, however, manifests itself as a divergence proportional to In^ , as can easily be verified 
using the results of Appendix B.1.3. This behavior is demonstrated in Fig. 4.10. 

The order g^ perturbative expression for the zero-temperature pressure, Eq. (|3.27|) . is ana- 
lyzed in Fig. 4.11, where it is plotted as a function of /i. Comparing its behavior with that 
of the dimensionally reduced 0{g'^) result in Fig 4.9, we see that as the temperature decreases 
the high-T expansion very smoothly approaches the zero-temperature one, until the unphysi- 
cal logarithms become dominant at small T. This indicates that our order g^\ng result may 

■^In Ref. j48j the breakdown of dimensional reduction has been studied more closely in the limit of large n/ 
with the result that its onset seems to take place already, when the ratio fx/T becomes of order 10. 
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Figure 4.11: Left: the scale dependence of the zero-temperature resuh of Eq. H3.27() . Right: the 
behavior of the order and g'^ perturbative results at T = and the order g^lng result at 
T = {).2Tc. The scale parameter is chosen as in Eq. ()4.4() . 

actually be applicable even at very low temperatures and high values of /x/T. Remembering 
on the other hand that the critical temperature for the meltdown of the color superconductor 
is of the order O.lTc, it would be extremely tempting to argue that the dimensionally reduced 
result is valid almost throughout the quark-gluon plasma phase. We must, however, exercise 
caution in drawing these conclusions until the order g'^ result applicable throughout the whole 
phase diagram is available. This computation obviously needs to be performed for the sake of 
completeness alone, even though it now seems that its numerical effect on our present results 
may well turn out to be of negligible magnitude in most parts of the (/i-T)-plane. 

To finish the chapter, let us still briefly dwell upon the properties of the zero-temperature 
expansion of the pressure. As we only have the first three perturbative orders available, it is 
somewhat difficult to quantitatively analyze the convergence properties of the series. A compar- 
ison of Fig. 4.11 with the corresponding ^ = plots of Figs. 4.1 and 4.2 however indicates that 
the situation at T = is considerably more favorable than at zero chemical potentials — with 
the possible exception of scale dependence. The magnitude of the 0{g'^) term in the expansion 
is nevertheless considerable, and the extension of the result to further perturbative orders is 
therefore certainly well-motivated. At the same time we must of course remember that in these 
calculations we are always developing perturbation theory around a physically incorrect ground 
state, which includes no Cooper pairs. 
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Chapter 5 

Conclusions and acknowledgements 



QCD is the quantum field theory that describes one of the four basic interactions of nature, 
the strong nuclear force. At extremely high temperatures and densities a strongly interacting 
system undergoes a phase transition from the hadronic phase into a new form of matter called 
the quark-gluon plasma, where the constituents of nucleons and other hadrons are liberated from 
their confinement and act almost freely like an ideal gas. In view of the on-going experimental 
efforts to detect and study this exotic phase of QCD it is essential that we continuously thrive 
to learn more about its properties. Analytic calculations in the theory should therefore be 
pushed forward, as they provide reliable and easily tractable results free from the practical and 
theoretical restrictions of numerical studies, such as lattice simulations. High-order perturbative 
computations, which provide a conceptually simple and efficient means to access the properties 
of quark-gluon plasma at very high temperatures and chemical potentials, certainly fall into this 
category. 

The work described in the present paper has dealt with the perturbative determination of the 
most fundamental quantity describing the equilibrium thermodynamics of quark-gluon plasma, 
the pressure of QCD at finite temperature and quark number density. After reviewing the 
necessary theoretical and mathematical background we proceeded to analyze the derivation of 
the most recent perturbative results for the quantity, and then moved on to compare them 
with up-to-date lattice data. As expected, our studies showed that the convergence properties 
of high-temperature perturbation theory improve considerably as we move away from the line 
fj, = 0. After analyzing the behavior of the different perturbative results for the pressure at 
small temperatures we then came to a somewhat surprising conclusion: the 0{g^ln g) expression 
derived for the quantity using three-dimensional effective theories in Ref. appears to be 
applicable almost throughout the quark-gluon plasma phase. At present this statement is, 
however, still of highly speculative nature at finite values of n/, and the existing order g'^ result 
for the zero-temperature pressure should therefore definitely be generalized to cover the region 
of small but finite temperatures as well. 

Taking into account the behavior of the pressure's high-temperature expansion at vanishing 
chemical potentials, it nevertheless seems to us that the most pressing challenge in the field is the 
determination of the 0{g^) correction to this quantity. At this order perturbative QCD, how- 
ever, runs into serious infrared problems, for which the only cure comes from non-perturbative 
methods, such as three-dimensional lattice simulations. As the results of four-dimensional lattice 
studies are on the other hand at present still a long way from covering the region of large T 
and our conclusion is that in the near future the most fruitful approach to high-temperature 
QCD may well be a combination of perturbative and non-perturbative techniques. First steps 



36 



in this direction have already been taken on both fronts [31 133 EDI- 
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Appendix A 

The matching coefficients 



Due to small notational differences between our treatment and that of Ref. f2| regarding the 
definition of the matching coefficients aui, 0(M2 and og, we include here a list of all coefficients 
required to obtain the pressure to 0{g^lng). In terms of the H functions introduced at the end 
of chapter 2 their explicit expressions read 
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+ 24 



52 H(3, z) + lUifL H(2, z) + (17 - 92/2^) H(l, z) + 4i/i H(0, , 



+ CfTp ( ^ (1 + (35 + 332/2^) - 24 (l - Afl^) 



144 



+ r| 



12i/i H(2, z) - 2 (1 + 8/i') K(l, z) - i/i (1 + 4/2^) H(0, 



^ (1 + 12/2^) (5 + 12/2^) In^ + 1 + 47 + 8 (7 + I27) fc^ + 112/2^ 



15 C(-3) 3 ^^ + ^^^ ) 



96 



8 H(3, z) + 12i/2 H(2, z) - 2 (1 + 2/2^) H(l, z) - i/2 H(0, , 



f fg 



+ 288T|— J^i 2(1 + 7) /2//2g- k3,z/ + Z3) + H(3,z/ + zp 



+ 4i/2/ (h(2, Zf + + H(2, Zf + z*)) - 4/2^ ^^(1, z/) - (/2/ + flg f H(l, z/ + 
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(fit - p-sf m,^i + 4)- ^vffii 2/) 
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3 Uf 



m A + C'(-i) 



AttT C(-1) 
r^l^ (1 + 12/^2) (2111^ +1) +24H(l,z)j|, 



? \Ca\ 22 In-— + 

CaTp ^2 (7 + 132fl^) ln|^ + 9 + 132/x2 + 87 + 4K(z)^ 
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^|cA(^221n 



e^A 
"4^ 



V 47rr 
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Appendix B 

Diagrammatic methods 



Up to this point our approach to perturbative QCD has been somewhat hmited, as the most 
comphcated computational aspect, the evaluation of multi-loop Feynman diagrams, has been 
left completely outside the treatment. In this appendix we will introduce to the reader the 
analytic techniques used in the articles ^ El IH] by presenting a few simple, but very detailed 
examples of our computations. In this process we will try to avoid a repetition of the material of 
the original papers, and rather concentrate on new and non-trivial aspects of the computations, 
such as the analyticity properties of the special functions encountered. 



B.l Three-loop diagrams in the full theory 

To illustrate the computational methods used in Ref. f2| in determining the strict perturbation 
expansion of the QCD pressure up to three-loop order, let us consider as an example the evalua- 
tion of the last diagram of Fig. 3.1 c, the fermionic 'Mercedes'. We begin the calculation by first 
reducing the diagram into a sum of master sum-integrals and then proceed to evaluate these 
functions. The computation will be performed in the Feynman gauge, where the value of the 
gauge parameter has been set to ^ = 1. 

B.1.1 Reduction to master sum-integrals 

Applying the finite-temperature Feynman rules of section 2 to the Mercedes graph, denoted here 
by Im, we straightforwardly obtain 



abc 



JjPQR} 

{P-Qf{Q-Rf{R~Pf 

where a sum over flavors has been suppressed and the sum-integration measure is as defined in 
Eq. (|2.84|) . The color sum is easily performed here by exploiting the antisymmetricity property 
of the structure coefficients f"'^'^, 

= '-cIaCaTf, (B.2) 



40 



and simplifying Eq. ()B.1|) further we eventually get 

Here we can at once take the trace over the gamma matrices by applying standard formulae 
(not specific to 4 dimensions) and get as the result a combination of metric tensors, which in 
our case are Kronecker (5-symbols. The subsequent tensorial contractions then become trivial 
and produce 

Ij^ = 48{l-e)dACATF9^^ P ■ (P - Q) Q ■ R 

^ ' T{PQR}P^Qm^{P-Qf{Q-Rf{R-Pf' ^ ' 

where the explicit appearance of e signals that we are working in 3 — 2e space dimensions. 
Combining the scalar products of the momenta into squares by using trivial identities such as 

P.Q = -i ((P-Q)2_p2_g2)^ 

we can then cancel two propagators in each term and finally get 

rO ^0^ - 1 



Im = 12(1- e) dACATpg^ (1^-1^)7 + - Mo,o , (B.6) 



where we have in analogy with Refs. j^J 1211 defined 

I^i - (B.7) 

^ " t{PQ} P^Q^ {P-Qf ^^'^^ 

Moo = T ^ o- (B.IO) 

nPQR} P^QHP - Rf {Q - Rf ^ ' 

Next we move on to consider the evaluation of these functions, which form an illustrative 
subset of the class of all one-, two- and three-loop master sum-integrals appearing in the covariant 
gauges of QCD. All sum-integrals required to determine the QCD pressure up to order g^\n.g 
have been evaluated in Ref. for /u = and in Ref. [2| for n 0. In the present work we will 
keep the chemical potentials arbitrary, and start our treatment from the fermionic one-loop case 
if. Our approach has been significantly motivated by the fj, = work of Ref. 0. 

B.1.2 Evaluation of the sum-integrals 

Written out using the integration measure chosen, the definition of if reads 



' - f{P} P' J (27r)3-2e p2 + ((2n + l)T - .^)2 • ^ • ^ 

After the trivial evaluation of the momentum integral (see e.g. the appendix of Ref. |B1]) we get 



= \^'lt:^ A^-r f: [{{2n+l).T-^^.fY'-\ (B.12) 

V^'^ ) n=—oo 



1/2-e 
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where the sum bears a close resemblance to the definition of the generalized Riemann ("-function, 

The sum-integral can therefore clearly be written in the form 

X (C(-1 + 2e, 1/2 - ifi) + C(-l + 2e, 1/2 + i/i)) , (B.14) 
which after an expansion in e — performed using the standard formulae of Ref. — gives 

ii = -^|l + 12A' + 2e (l + 12/i2) (^l + ln^^ +12H(l,z) |. (B.15) 



Here we have switched to using the MS scale parameter A and now encounter the H-functions 
for the first time. The bosonic sum-integral if is evaluated in an analogous manner. 

As we move on to consider the two-fold sum-integral T — first solved in Ref. j2,Sj — the com- 
putational methods required get somewhat more complicated. The initial step in the calculation 
is to notice that the definition of the sum- integral, Eq. I|B.9|) . can be written in a simpler form, 
if we shift the integration variables by 

P P' = P-Q, (B.16) 

and define a 'polarization' function Ilf as 

We then get for the sum-integral 

r ^ t ^^^^i . = i^MP), (B.18) 



P{Q} P^Q^ {P + Qf % P^ 

and begin its evaluation by transforming the function Ilf into coordinate space at e = 0. 

As can be easily verified, a Fourier-transform in three dimensions gives for the fermionic 
propagator 



j3 iq-r -(|2n+l|7rT-i/^ sign(2n+l))r 

^ . (B.19) 



(27r)3 q2 + ((2n + 1) vrT - i^if 47rr 
Inserting the inverse transform into the definition of Ilf we then obtain (po = 27rmT) 

X exp[ — (|2n -|- 1\t:T — i/isign(2?7- -|- 1)) ri] 

X exp[- (|2n + 2m-M|7rr-i/isign(2n + 2m-M))r2], (B.20) 
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where the momentum integral can be easily performed and produces a Dirac (5-function. This, 
on the other hand, helps us to remove one of the r-integrals and gives 



nf(P) 



1 ^ r- 

—-^ / dVe*P-'exp[-(|2n+l|7rr-i/xsign(2n + l))r] 

^ ' n=-oo'' 

X exp[ - (|2n + 2m + lIvrT - i/xsign(2n + 2m + 1)) r], (B.21) 
in which the infinite sum over n is quickly performed with the help of the standard formula 



m 

n=no 



csch(a/2) ^^_„(„g_i/2) 



-a(ni+l/2) 



(B.22) 



Our result for the polarization function in three dimensions then finally reads 



nf(p) 



T 



(47r)^ J r 
with the dimensionless variables po and f defined by 



d r — ^ (|po| + cos{2fir) cschr) e 



-\po\r 



PO 



2ttT ' 
f = 2'KTr. 



(B.23) 



(B.24) 
(B.25) 



The insertion of Eq. (|B.23|) into the definition of T provides us a straightforward way to 
compute the sum-integral. Before doing this we however must find a way to regulate its obvious 
UV divergence, which otherwise would render the value of the graph infinite. The most straight- 
forward way to proceed is to analyze the behavior of Ilf at large external momenta and then 
subtract the leading order terms from the sum-integral. The subtracted part namely becomes 
finite and can be evaluated at d = 3, while the rest is easily dealt with while keeping e nonzero. 

The large-P behavior of the polarization function is most easily obtained through the standard 
procedure of converting the go sum in Eq. ()B.17|) into a contour integral (see e.g. Refs. 
using the formula 



°° rioo 
^„ J —ioo 



dpofiPo) + f^(ipof{Po) 



too 

-ioo+fi+e 



(B.26) 



ioo+fi—e 



dpon{fi-po) }fipo 



where we have denoted 



n{x) 



= (2n + l)7rr - i/i, 
1 



(B.27) 
(B.28) 



and the integration contour C is defined by 

ioo 



C = 







-too 



Po 



(B.29) 
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Concentrating first on the divergent vacuum (T = /j, = 0) part of Ilf, represented by the first 
term in Eq. HB.26|) . we straightforwardly obtain 



n(o)(p) = I 



d^-^'q 1 



K^' I dx 



1 f Ai-2e 







d^-^'q 1 



(27r)4-2^ (g2 + xP2 + 2x q • P)2 



A2T(e) , 1 
dx- 



(47r)2-^ 7o (x(l-x)p2). 

- 1 r(.)r^(i-6) /A^Y _ /T^Y 

where /Jq clearly contains an 1 /e divergence. We have suppressed the subscript f, as the vacuum 
part of the polarization function is independent of the bosonic or fermionic nature of Q. 
Let us then move on to consider the matter part of ITf and define 

f{m) = ^ 9 0- (B-31) 

(i(<7o+Po))2-(q + p)' (^'Zo)2-q2 



Following the notation of Eq. (|B.26|) . we now get 

iif\p) = nf(p) -n(°)(p) 



dgo/(go) 



{90} 

^2e r- ^3-2e 



2m J (27r)3-2^ y^j^ 
- \ / dqon{qo-fi)+ dqo n{fi - qo) i f (qo) ) , (B.32) 

in which the go-integrals are trivial to perform using the Residue theorem. After some straight- 
forward simplifications this produces 

Uf\p) = A'^ J 1 (^(/x -q)+ n{q + + s\gn{q - ,,)n{\q - ^,\) 

1 1 

+ 



^2. 



^Z-2e^ 1 



(27r)3-2e 2q 



n{q + fj,) + n{q - fi)^ 



. {q - ipoY - (q + p)^ ^ {q + mY - (q + p)0 ' (b.33) 

where we have used the trivial identities 

niq-ipo) = n{q), (B.34) 
niq-fi) = 6'(/i-g)+sign(g-/i)n(|g-/i|). (B.35) 

(T) 

The leading large-P behavior of is now available through a simple power series expansion 
of the integrand of Eq. HB.33|) , as the functions n ensure that only the region of small q contributes 
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to the integral. This on the other hand helps us to immediately take care of the angular integrals 
and leads to the result 



nf)(p) 



A- 



2e 



X 



(47r)3/2-.r(3/2 - e 
f 2 



dqq 



l-2e 







n{q + ii) + n{q - ix)^ 



+ 



' 3-2e 



q^ + O 



p6 



where the remaining g-integrals are expressible in terms of the polylogarithm function 



OO I, 



E 

k=l 



r(^) 



.u-1 



dt- 



(B.36) 



(B.37) 



Our final result for the leading large-P behavior of the matter part of the polarization function 
then reads 



with the coefficients /3 defined by 



z vr 



-3/2+e 



r(3 - 2e) 
r(3/2-e 



(Li4-2e(-e'^/^)+Li4„2e(-e- 



(B.38) 

(B.39) 
(B.40) 



As can be easily verified starting from the definition of Ilf, we could have immediately written 
down the first coefficient in the form /3i = 2X^^ which actually is enough for our present purposes. 
We, however, wanted to present here the general scheme used in analyzing the large-P behavior 
of different polarization functions and therefore extended the treatment to the next order. The 
higher order terms are frequently needed at three-loop level. 

Returning back to the computation of T, we now see that the explicit subtraction of the 
leading order terms of Ilf from Eq. ()B.18|) is enough to render the sum-integral finite. We, 
however, still need a coordinate space representation for the regularized polarization function at 
e = 0, which is straightforwardly obtained and reads 



nf(P)-n(o)(p)-nfuV(^) 



(B.41) 



(4vr) 



dV- 



cos(2/xr)csch r — - + [-+ 2^ ) r 



1 



V360 3 3 



-\po\r 



In deriving this result, we have employed e.g. the identities 



Li4( 



Li2(-e^') + Li2(-e-" 



vr 



+ Li4( 



vr 



360 12 



24 



(B.42) 
(B.43) 



which are special cases of a more general relation between the polylogarithm function and the 
Bernoulli polynomials (see Ref. [HSI)- As a curiosity, it should be noted here that despite the pole 
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in the coefficient /3o, even the Fourier transform of II'^'^^ is finite at d = 3 due to the appearance 
of an additional factor of l/r(e) in the transformation formula of l/{P'^y. 

We are now finally ready to evaluate the initial sum-integral. Plugging the above results into 
Eq. (|Bl8l) . we get 



T 



+ 



1 

T 



/r2 



(Any 



p2 

cos(2/ir)csch r — - + - + 2^' 
r V 6 



\PO\r 



+ 0(6), (B.44) 



where the two first terms are easily computable by methods analogous to those used when 
considering I^. For the finite third term the evaluation of the momentum space integral and 
the subsequent po sum on the other hand produce 



T 



(An)- 
(4^ 



1 



d^. 



-2|pok' 



Po 

°° , cothr 

dr 

r 



cos(2^r)cschr 
cos(2/if)cschf - 



1 



+ 



+ 



,-bol'- 



+ 2/2^ 



cos(2^r)cschr 



+ 



(B.45) 



where the final coordinate space integral is UV convergent but has an IR divergence due to 
the zero mode in the bosonic frequency sum. We could easily remove this by separating the 
zeroth Matsubara mode from the sum, but need not actually do that; the singular part namely 
originates from a term proportional to Jg°° dr r^, which automatically vanishes under dimensional 
regularization. We can therefore instantly proceed to evaluate the integral, which is most easily 
accomplished by regulating each term of the integrand separately with a factor r". In the end 
we will sum the different pieces together and proceed to the limit a — > 0+. 

The last two terms of Eq. (|B.45|) are easily dealt with by applying the relation 



dx x^coth X 



2-^r(i + z)c(i + z), 



(B.46) 



while with the first one we first need to perform a partial integration and then use the straight- 
forwardly derivable identity 



dxx^e^^^cschx = 2-^r(l + z)C(l + 1/2 - i/5/2). 



(B.47) 



When adding the three parts of Eq. HB.45|) together, all terms singular in a cancel and we are 
left with an expression, where we may easily proceed to the limit a = 0. The inclusion of the 
first two terms of Eq. (|B.44|) to the sum then finally leads us to the result 

T^p,^ /I . _ A \ 2i 



T 



1 ^ , A 



H(0,z) 



(B.48) 



for the entire sum-integral. 

For three-fold sum-integrals the calculations become still somewhat more evolved [21121) but 
the general features of the above coordinate space method still apply; in particular, we always 
end up evaluating one-dimensional hyperbolic integrals, which lead to results containing different 
combinations of the J^-functions. We will therefore not present a detailed account of the com- 
putation of the sum-integral Mo,o here, but merely refer the reader to Ref. El. In the following 
we instead concentrate on analyzing the properties of the special functions encountered. 
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B.1.3 Properties of the K-functions 

As can be seen from the explicit form of the matching coefficients listed in the appendix A, the 
only special functions that appear in the perturbative expansion of the QCD pressure up to 
order g^ln g are the ^^-functions. These were defined in chapter 2 by 

H(n,u') = C'i-n,w) + {-l)''+'^C'{-n,w*), (B.49) 
i<{w) = "^{w) + '^{w*), (B.50) 

where n G {0, 1, 2, 3} and the arguments of the functions read either 

w = z = 1/2 — z/i or (B.51) 

W = Zf + Zg = 1 - i{flf + fig). (B.52) 

We will now examine the behavior of these functions and in particular show how one can derive 
generalized power series expansions for them in the limits of small and large fi. These results have 
played a crucial role in the derivation of analytic expressions for the quark number susceptibilities 
in Ref. |^ and for the T = pressure in Ref. [2|- 

Our first task is to find convenient (integral) representations for the C- and ^'-functions. The 
generalized Riemann ("-function is usually defined by the infinite sum 

OO ^ 

C{x,w) = y2-—— (B.53) 



or the integral j66j 



1 I p-->"t 
Cix,w) = ^ / dtt^-'-^-^. (B.54) 



Using the latter it is especially straightforward to analytically continue^ the function to the 
whole j;-plane excluding the point x = 1, if we assume the parameter w to satisfy the condition 
Re(t(;)> 0. For the digamma function ^ we will use a similar representation |66j 



= -7+ / dt— —, (B.55) 

Jo 1 - e * 

where Re(t(;)> is again implicitly assumed. 
Small p, 

Let us begin our analysis from the limit of small fl, where we must expand the integrals of 
Eqs. (|B.54|) and (|B.55|) in powers of lm{w). For the (" function this is easily achieved via an 
expansion of the denominator of the integrand in powers of e~*, which assumes that the factor 
^x-i regulates the divergence at small t. If we denote Re(i(;) = wq, Im(^i;) = wi, we get 

OO 



-1 ^ noo 



n=0 

^ CXD OO / , \^ /'OO 



T(x) k\ 

^ ' n=Ofc=0 



^Recalling the analytic continuation of the V function onto C\(— N) using its integral representation. 
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= f;(^^!"^)c(A; + x,«;o)(-iu;l)^ (B.57) 
fc=o ^ ^ 

where the binomial coefficient in the final expression reveals that the result could as well have 
been derived starting from the sum representation of A derivative of the above expression 
evaluated at negative integer values of x then gives us the desired results for tt;). Just to 
demonstrate this, for n = 1 the first few coefficients of the expansion read 

^^(1,^) = -— (ln2-^|^) -(l-21n2-7)/22 + 0(//^). (B.58) 



12 V C(-l) 

For the functions 'ii{w) and ^ the calculation becomes somewhat more complicated. Intro- 
ducing an IR regulator to Eq. ()B.55|) and expanding the denominator as before, the integrals 
corresponding to the two terms of the integrand can, however, be analytically computed and 
produce a finite result when added together. We are then able to expand also the function ^(z) 
to an arbitrary order in fi, and obtain 

= -27-41n2 + 14C(3)/22 + C)(^4^_ ^g^g^ 



Large fi 

Moving on to the limit of large ft, we want to expand the H-functions in powers of or 
equivalently the integral representations of Eqs. ()B.54|) and (|B.55|1 in powers of l/lm(w). This 
is most easily accomplished by noticing that a large value of p, causes the dominant contribution 
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to the integrals to originate from the region of small t, which means we can expand the integrands 
in powers of t with the exception of the factor e~™*. We thereby get 

({x,w) = —— / dtt"-ie-™*(t-i + l/2 + t/12 + 0(t2)) 

r(x) Jo 

X — 1 \ 2iw 12[iw)'^ ) 
which after some trivial algebra leads to 

= -^2(^ln/i-l)-l(ln/z + l)+0(^i^) (B.61) 

for z\ and similarly for the other H(n, z)-functions. 

Analogously to the limit of small /i, the determination of the behavior of the functions "^{w) 
and ^ at large \i requires the introduction of an additional factor t° to the integrand of Eq. HB.55|) 
to regulate the divergent behavior of the two terms at t = 0. In other respects the computation 
however proceeds in a manner highly analogous to that presented above and in the end gives us 



= 21n/i-li^ (B.62) 



The general behavior of the functions z) and is depicted in Figs. B.l. and B.2. 



B.2 Four- loop diagrams in the effective theories 

In Ref. |3I it is described how all four-loop vacuum diagrams appearing in three-dimensional 
massive field theories can be numerically computed to a high order in e = (3 — d)/2. For 
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Figure B.3: Diagrams evaluated in the present paper. 



these purposes we have in addition derived previously unknown analytic results for a small 
class of four-loop master diagrams, whose structure is simple enough to allow for their analytic 
evaluation. As this analytic work is only briefly quoted in Ref. S:, we will here explore the 
subject in more detail and introduce the methods used in the computations. An even more 
detailed treatment of some of this material can be found in Ref. |67j . 

The diagrams, whose analytic evaluation turns out to be possible, fall into two general cate- 
gories: graphs with a small number of vertices can be evaluated by coordinate space methods 
(see e.g. Refs. |681 169j ). while the ones containing only simple one-loop sub-diagrams are most 
easily dealt with directly in momentum space. Our treatment is therefore too divided into two 
parts, of which the first deals with coordinate space computations and the latter with momentum 
space ones. The diagrams to be covered are shown in Fig. B.3, where also the numbering of the 
propagator masses can be found. We will keep the values of these masses arbitrary throughout 
our treatment, which enables us to afterwards easily obtain results for diagrams with some prop- 
agators raised to higher powers through differentiation. The momentum integration measure we 
will be using in our calculations is the MS one defined at the end of chapter 2. 



B.2.1 Coordinate space calculations 

When the number of vertices in a diagram is small, it often proves useful to convert the momen- 
tum space integrals into coordinate space ones by taking Fourier transforms of the propagators. 
The success of this method is based on the obvious fact that the number of coordinate space 
integrations needed in the end will always be one less than the number of vertices in the graph. 
The coordinate representation of the scalar propagator is, however, fairly complicated in an 
arbitrary dimension, and angular integrals furthermore often prove to be surprisingly difficult 
to perform. This makes the application of the method increasingly difficult when the number of 
vertices arises. 

The basic setup 

The most fundamental quantity required in coordinate space calculations is the Fourier transform 
of a massive propagator, for which a straightforward computation produces 

D(x,m) ^ A2^ / 4^ ^— ^ = - Ki/2_,(mx), (B.63) 

J {2tt) p'^ + m^ (27r) ' ^ ^ x J ' 

with K denoting a modified Bessel function. An expansion of this result in powers of x then 
gives 



^ ' ' 47ri-T(l/2) \ 2{l + 2e) 

- A=^T (-1/2 + e) (i + oi {mxf ) ) (B.64) 

(47r)i~^r(-l/2) V ^ ^ > 



50 



l-2t 



= ~" a [mx) 
(3 (l + o({mxf 



-i+2e I ^ ^ yrnx) 



2(l + 2e) 



+ O ( {mx 



where the coefficients a and (3 read 



a 



A^T(l/2-e) 
47ri-^r(l/2) ' 
A2T(-l/2 + e) 



(4vr)i-r(-l/2) 
Correspondingly, an e-expansion of D produces 



D (x, m) 



47rx 



1 + e( In^^!^ + e2™-'Ei (-2mx) ^ + O (e^) 



where Ei denotes the exponential integral function 



Ei(2;) 



(B.65) 

(B.66) 
(B.67) 

(B.68) 

(B.69) 



and the relations j66j 



Kl/2 {x) 




dK^ (x) 
du 



v=l/2 



2x 
vr \i/2 
2i 



Ei (-2x) 



(B.70) 
(B.71) 



have been used. These results will be frequently used in the following. 
Diagram d 

The basketball-type diagrams containing only two vertices are particularly well suited for inte- 
gration in coordinate space due to their trivial angular structure. Writing the propagators in 
terms of their Fourier transforms we easily obtain for the graph d, 



Id 



A 



-lOe 



A 



2e 



(47re-'>) 

A-i 
(47re-T) 



(47re-'>) 

4 



i=l 



27r3/2- 



r(3/2-e) 
2^ 27r3/2-^ 



dx X 



2-2e 



x,mj) 



r(3/2 



(B.72) 



where the integral X clearly contains a singularity at the origin signalling of the ultraviolet 
divergence of the diagram. In order to separate this divergent part we introduce a small positive 
parameter r to divide the radial integration path into two subintervals [0, r] and [r, oo[ while 
keeping in mind that in the end we will sum the two contributions together and proceed to the 
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limit r ^ 0. On the first interval one may use the small-x expansion of the propagators, whereas 
on the latter we can simply set e = 0. 

Starting from the first interval we get by applying Eq. ()B.65|1 . 

5 



Ii = dx x^~^'TTD(x,mi) 
•^0 i=i 

■^0 i=l 

-I 



a [miX) 
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where the final integrations are trivial to perform and lead to 
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On the second interval the calculation is even more straightforward, as at (i = 3 we only need 
to evaluate the integral 



l2 



1 



(4^) 



dx 



with M defined by 



M 



E 

i=l 



(B.75) 



(B.76) 



After several partial integrations and the dropping of terms that vanish in the limit r ^ we 
eventually obtain 

-Mr ^e-Mr M^e-^'^Hu (Mr) 



To 



1 
1 



2r 



J_ _M_ 1 / 3 _ 
(47r)^ V2r2 ~ T 2 V2 ~ ^ ' 
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When adding up Xi and X2, all terms divergent in the limit r ^ cancel. An expansion in e 
then gives for the whole integral the simple result 



Id 



2(47r)4 [8e 
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A 
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3 ^ A 
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(B.78) 
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X 1 y 






Figure B.4: The structure of the diagram e in coordinate space. 

where we have finahy proceeded to the hmit r = 0. When all masses are equal, = mVi, we 
in particular get 

/,(m,, = m) = ^ fl+SlnA+iS-f In2-^ln5y (B.79) 



Diagram e 

The coordinate space evaluation of diagrams containing more than two vertices is generally very 
difficult to approach, but may anyhow prove successful in certain special cases. One such case 
is the four- loop graph e, which after a Fourier-transform reads 



h = (77^— ^1 / d3-2^x / d3-2^y D(|x-y|,mi 



A- 

(47re-T) 

X D{x,m2)B{x,m3)D{y,m4,)B{y,m5)D{y,me) , (B.80) 

as portrayed in Fig. B.4. Here we can at once perform the angular integrations by applying the 
identity [HHl 



/ 



d^P^¥^ = {2nr''-^^, (B.81) 



where A = [d — 2)/2 and r > p. This on the other hand gives us 



e - 17^^3^J dx dy x^~^'y^-^'l[Bix,mi)llDiy,„.,, 



L = I J- I / dx / dy x"""V-"M I Dfx,mi) M Df?y,m, 

(47re 



i=2 i=4 

X (D (x, mi) G (miy) 6* (x — y) + D (y, mi) G (mix) 9 {y — x)) 
in which we have defined the function G by 

25/2-e 3-2.^^-4. J U) 

For future purposes we note that at small values of its argument this function has the expansion 
G(x) = ^^^2i} + 0{x^)) - + (B.84) 



and with e set to zero it reads 



G(x) = (47r)2^i^^. (B.85) 
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Figure B.5: The division of the xy-plane into 6 subregions in the evaluation of the integral of 
Eq. (lEHS). 



Let us then move on to evaluate the integral I. In analogy with the computation of the 
basketball diagram we again introduce a small positive parameter r to divide the integration 
plane into different regions in order to isolate the singularities of the graph. This division is 
depicted in Fig. B.5, and it is obvious that at least the two regions E and F will be safe from 
UV divergences meaning that we can there set e = 0. In the end we will again sum the different 
contributions together and proceed to the limit r — > 0. 

In the computations to follow we will frequently apply the notation 

rrii + rrij + + ... = Mij^^^^, (B.86) 
X xm'' X ... = M*^*^-. (B.87) 



Regions A and B 

In the regions A and B the coordinates x and y are confined near the origin as x,y < r. Denoting 
y = tx and using the small- a; expansions of the different functions we easily obtain 



Ia = [ dx r dy x^-^'y^-^'G{miy)f\D{x,mi)f\D{y,m^) 

Jo Jo .^4 

= Ka^ dx dt + 0(r) = Ka^^+0{r), (B.88) 

Jo Jo 32e 



and similarly 



Tb = dy dx x2-V~''G(mix)D(y,mi)rrD(rE,mi)TTD(y,mi) 
As we will in the end take the limit of r = 0, this accuracy suffices to us at present. 



Region C 

In the region C — which turns out to be the most problematic one — the integration variable 
y in confined near the origin, while x on the other hand ranges from r ro oo. Proceeding in the 
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same fashion as before, we get 



poo rr ^ ^ 

Ic = / dx / dy x2-V"''G(miy)nD(x,mi)J]D(y,m,) 

i=l i=4 

f'oo ^ pr 

= na^ dxx2-2^JjD(x,mi) / dy y-^+^' + O {r) 
Jr -^^ Jo 

= Ka^ — / dx x^-'^^WV)(x,mi) + Oir), (B.90) 
4e Jr 



1=1 

where the remaining x-integral must be computed to 0{e) in order to obtain also the finite part 
of Tc- This can on the other hand be achieved through an apphcation of the e-expansion of D, 
Eq. (|Rfi8|) . which gives us 

/ dx x^-^' JJ D (x, mi) (B.91) 
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(47r) 



2\3 3 



i=l 

but leaves a number of non-trivial one-dimensional integrals to be computed. Performing several 
partial integrations all of these can nevertheless be evaluated analytically, and we finally obtain 

+ \ (In (Mmr))" + ^ (In (2mir))" + 7I11 (2mir) + Lij (l - ^ ) ] | + 0(r) 

i=l ^ ^ / ) 

as the contribution of the region C to the integral I. 
Region D 



roQ 3 6 

Id = dx dy x^y^G(mix)D(y,mi) JjD(x,mi) JjD(y,^ 

^0 Jr • r, ■_, 



It can be easily verified that the region D produces no divergences, so we may set e = in 
the beginning. Then all the required integrals reduce to familiar types and we easily obtain, 
following the above treatment, 

3 6 

i=2 j=4 

= 2 / dx sinh (mix) / dy 7: 

(47r) mi Jo X Jr r 

= +0(r). (B.93) 

{Any 

Regions E and F 

The regions E and F contain no points near the origin and will therefore also give only finite 
contributions to the integral X. Proceeding as before we get for E, 

/oo poo 3 6 

dy / dx x^y^G (miy) J] D (x, m,) J] D (y, m,) (B.94) 
•^y i=i i=A 
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and similarly for F, 
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Here the remaining integrals are straightforward to evaluate by employing partial integrations 
and in the end produce the results 
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All that is left to do is to add up the contributions of the different regions to the initial integral 
and proceed to the limit r ^ 0+. This finally gives us as the result for the whole diagram 
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Special mass configurations 



Important special mass configurations include the identical mass case, m^ = m Vi, as well as the 
one with m3 = 0, mi = m2 = m4 = = uiq = m. Using the easily verifiable dilogarithmic 
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relations 
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(B.102) 

where especially the remarkably simple form of the latter result is worthy of some attention. 



B.2.2 Momentum space calculations 

Diagrams containing arbitrarily many vertices but only trivial one-loop subdiagrams are most 
easily handled directly in momentum space. The number of such graphs is, however, small and 
analytic results are so far available only in a few special cases. We will now briefly address the 
computation of three four- loop diagrams, for which at least simple integral representations can 
be derived using this method. 



Diagram a 

The triangle diagram has previously been studied both in three and four dimensions |7n| I71j 
with the result that its divergent part has been successfully obtained. Here our objective is, 
however, to compute its finite term as well, and we start this by defining a scalar two-point 
function 



fl{p,mi,m2) = [ j 3 ^ 



(B.103) 



analogous to the polarization function of appendix B.l. This allows us to write the diagram in 
the form 



la 



n(p, mi, m2)n(p, ms, mi)Il{p, m^, mg). 



(B.104) 



which serves as a useful starting point for our analysis. 

The momentum space integral appearing in the definition of 11 can be easily performed by 
introducing a Feynman parameter x to combine the two propagators in the integrand (see also 
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Eq. dOfl V We thereby get 
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where the function B possesses the easily verifiable properties 
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which after a partial integration in u gives 
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Here the 1/e part of the integral has been separated, which means we can expand the factor 
in e. Using Eqs. (|B. 106(1 and l|B.107() this finally leads us to the result 
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where terms of order e have been dropped. 

The final result for the triangle diagram is now available through a simple e-expansion of 
Eq. (|B. 110(1 . which gives us 
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Here the remaining integral can be analytically performed only in certain special cases, which 
fortunately include the one with identical propagator masses (mj = m Vi), 
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Diagrams b and c 

The graphs b and c are clearly finite in three dimensions. Using the above results for the function 
n we may almost instantly write for them the integral representations 
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X arctan ( ) arctan ( ) arctan ( ) , (B.113) 
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(B.114) 



The analytic evaluation of these integrals has, however, not been possible even in the identical 
mass case, so at present these expressions only serve as sources for accurate numerical results. 
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